International Journal of Modern Physics D, 
© World Scientific Publishing Company 



THE i?-MODE INSTABILITY IN ROTATING NEUTRON STARS 



NILS ANDERSSON 

Department of Mathematics, University of Southampton, Southampton S017 IBJ, United Kingdom 

KOSTAS D. KOKKOTAS 

Department of Physics, Aristotle University of Thessaloniki, Thessaloniki 54 006, Greece 

Received (received date) 
Revised (revised date) 

In this review we summarize the current understanding of the gravitational-wave driven 
instability associated with the so-called r-modes in rotating neutron stars. We discuss 
the nature of the r-modes, the detailed mechanics of the instability and its potential 
astrophysical significance. In particular we discuss results regarding the spin-evolution 
of nascent neutron stars, the detectability of r-mode gravitational waves and mechanisms 
limiting the spin-rate of accreting neutron stars in binary systems. 

1. Introduction 

The first decade of this new millcnium holds great promise for gravitational 
physics. The hope is that forty years of development will come to fruition as the 
new generation of gravitational-wave interferometers (LIGO, VIRGD, GEO600 
and TAMA300) come online and reach their projected sensitivitieslll. These de- 
tectors should (finally!) open a new window to the universe and make the long- 
heralded field of "gravitational-wave astronomy" a reality. For gravitational- 
wave theorists this means that decades of modelling will finally be tested by ob- 
servations. The construction of the new detectors is accompanied by attempts to 
model all likely gravitational-wave sources in appropriate detail. This is highly 
relevant since reliable theoretical templates, against which one can match the 
noisy data-stream of the detectors, are needed if we want to extract reliable 
astrophysical information from future gravitational-wave data. 

In this article we survey the territory where stellar pulsation theory meets 
gravitational-wave astrophysics. There has been a resurgence of interest in this 
research area (that dates back to seminal work by Thorne and his colleagues 
in the late 1960dj'BQ) and our aim here is to provide an overview of the recent 
developments. We will focus our discussion on the gravitational-wave driven 
instability of the so-called r-modes. Since its serendipitous discovery a couple 
of years agoill, the r-mode instability has attracted considerable attention. In 



1 



2 The r-mode instability in rotating neutron stars 



this review article we aim to decribe most of the recent ideas and suggestions, be 
they speculative or not, and put them in the appropriate context. We will try to 
introduce the relevant concepts and results from basic principles, in such a way 
that the article should be accessible to readers with little previous knowledge of 
this field. We feel that this is an important task since it may help us identify 
the many outstanding issues that must be addressed by future work. 

The theory of stellar pulsation is richly endowed with interesting phenom- 
ena, and ever improving observations suggest that most stars exhibit compli- 
cated modes of oscillation. Thus it is natural to try to match theoretical models 
to observed data in order to extract information about the dynamics of dis- 
tant stars. This interplay between observations and stellar pulsation theory is 
known as asteroseismologyQ. It is interesting to speculate that the advent of 
gravitational-wave astronomy will permit similar studies for neutron stars. One 
would expect most neutron stars to oscillate during various phases of their life, 
and actual observations would provide invaluable information about the stars 



pactness, with a mass above 1.3Mq concentrated inside a radius of a mere 10 km, 
oscillating neutron stars could be interesting astrophysical sources of gravita- 
tional waves. But it is not clear that astrophysical mechanisms can excite the 
various oscillation modes to a detectable level. It seems likely that only the 
most violent processes, such as the actual formation of a neutron star following 
a supernova or a dramatic starquake following, for example, an internal nhase- 
transition, will be of relevance. The strengthening evidence for magnetarsau, in 
which a starquake could release large amounts of energy, is also very interest- 
ing in this respect. One can estimate that these events must take place in our 
immediate neighbourhood (the Milky Way or the Local Group) in order to be 
observable. 

This does not, however, rule out the possibility that neutron star pulsation 
may be detected by the new generation of gravitational-wave detectors. First 
of all, one should not discard the possibility that unique events in the life of 
a neutron star may excite the various modes to a relevant level. More likely 
to be observationally important is the long recognized possibility that various 
pulsation modes of rapidly rota.tin g; neu tron stars may be unstable due to the 
emission of gravitational radiationli3£3o. Should such an instability operate 
in a young neutron star it may lead to the emission of copious amounts of 



so-called r-modes, for which the gravitational waves have been estimated to 
be detectable for sources in the Virgo cluster (at 15-20 Mpc). If we suppose 
that most newly born neutron stars pass through a phase where this kind of 
instability is active, several such events should be observed per year once the 
advanced interferometers come into operation. This is a very exciting prospect, 
indeed. 
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Of particular recent interest is the instability of the 
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2. Nonradial oscillations of rotating stars 

A neutron star has a large number of families of pulsation modes with more 
or less distinct character. For the simplest stellar models, the relevant modes 
are high frequency pressure p-modes and the low frequency gravity g-modeg3. 
For a typical nonrotating neutron star model the fundamental p-mode, whose 
eigenfunction has no nodes in the star, is usually referred to as the /-mode. The 
/-mode has frequency in the range 2-4 kHz, while the first overtone lies above 
4 kHz. The 5-modes depend sensitively on the internal composition and tem- 
perature distribution, but they typically have frequencies of a few hundred Hz. 
The standard mode-classification dates back to the seminal work of Cowlingl3, 
and is based on identifying the main restoring force that influences the fluid 
motion. As the stellar model is made more detailed and further restoring forces 
are included new families of modes come into play. For example, a neutron star 
model with a sizeable solid crust separating a thin ocean from a central fluid 
region will have g-modes associated with both the coieand the ocean as well as 
modes associated with shearing motion in the crusttZHl3. Of particular interest 
to relativists isJJiLLCxistence of a class of modes uniquely associated with the 
spacetime itseliEjSEil; the so-called w-modes (for gravitational wave). These 
modes essentially arise because the curvature of spacetime that is generated by 
the background density distribution can temporarily trap impinging gravita- 
tional waves. The w-modes typically have high frequencies (above 7 kHz) and 
damp out in a fraction of a millisecond. It is not yet clear whether one should 
expect these modes to be excited to an appreciable level during (say) a grav- 
itational collapse following a supernova. One might argue that they provide 
a natural channel for the release of any initial deformation of the spacetime, 
but there are as yet no solid evidence indicating a significant level of if-mode 
excitation in a realistic scenarioil@0. 



2.1. Linearised equations of motion 

In this article we will only touch briefly on "realistic" neutron star models. 
Most of our discussion will concern (adiabatic) linear perturbations of relatively 
simple perfect fluid models. We will typically assume that the unperturbed 
pressure p and density p are related by a one-parameter equation of state, p = 
p{p). A useful way of writing this is to introduce an effective polytropic index 

by 

dlogp 



dlogp 



(1) 



In the special case of a polytrope F is, of course, constant and often expressed 
as F = 1 -I- 1/n. 

We want to consider a star rotating uniformly with frequency Q. In the 
frame of reference that co-rotates with the star the equations that govern the 
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fluid motion (the Euler equations) can be written 

dtu + u ■ Wu + 2n X u + n X (d X f) = -ivp- V$ (2) 

P 

where u is the fluid velocity and <I> represents the gravitational potential. Here 
it is worth noticing that the centrifugal term can be rewritten as 

nx {nxr) = -^V(l1 X f)2 (3) 

In other words, hydrostatic equilibrium corresponds to 
1 



-Vp = -V 
P 



= -V* (4) 



From this we can deduce that the level surfaces 

* = $ - sin^ l^l = constant (5) 

coincide with the isobars (and isopycnic surfaces). Consequently, once we have 
found a convenient way to describe the level surfaces of the "effective gravita- 
tional potential" 4" we have a complete description of our rotating fluid. In the 
particular case of slow rotation, the level surfaces can be described by introduc- 
ing a new variable a corresponding to the mean distance of a given constant 
^ surface to the centre of mass. This new variable is related to the standard 
spherical coordinates (r and 0) through 

r^a[l + e{9)]. (6) 

where a ranges from to i? (the radius of the corresponding nonrotating star). 
The advantage of this representation is that the pressure and the density can 
now be thought of as functions of a only. For a uniform density star (which 
incidentally is not a bad approximation for a neutron star) we have 

where P2{cos9) is the standard Legendre polynomial, and for polytropes one 
can calculate the required function e{6) from the results of Chandrasekhar and 
LebovitzEa. 

Having prescribed a rotating stellar model, we want to study small nonradial 
perturbations away from equilibrium. We will consequently assume that the 
perturbation is characterized by a displacement vector ^, which is small in a 
suitable sense, and yields the fluid velocity as 

Su = dt^ (8) 
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Given that the background model is stationary and symmetric with respect to 
the rotation axis, we can always decompose a pertm'bation into modes that 
depend on the azimuthal angle as exp(im(p) with m an integer. If we also 
assume that these modes have a harmonic time-dependence we can write 

where ujr is the oscillation frequency measured in the rotating frame. By recall- 
ing that the inertial frame is related to the rotating frame by 

d/dt^ d/dt + u-V ^ dt + nO^ (10) 

we see that the inertial frame frequency is given by 

LUi = LUr — mQ . (11) 

With the Ansatz (^) the linearised Euler equations are 

dtSu + 2nxSu^ -t^rC + 2it^rf^ X ^Vp - -VSp - V6<i> (12) 

P P 

where dp and Sp represent the Eulerian perturbations in the density and pres- 
sure, respectively, and the variation in the gravitational potential is Note 
that the centrifugal force only enters this equation through its effect on the 
background configuration. That this should be the case is easily understood 
from (U): Since we are working in an Eulerian framework the position vector r 
is treated as a constant. 

To fully describe the perturbations we also need the equation describing the 
conservation of mass 

dtSp + V ■ (pSu) =0 (13) 

or in integrated form; 

Sp + V-{pO^0 (14) 
as well as the perturbed Poisson equation for the gravitational field 

V2j$ = AirGdp . (15) 

We must also prescribe an "equation of state" for the perturbations, i.e. a rela- 
tion between the (Lagrangian) perturbations of pressure and density. Assuming 
adiabatic perturbations this relation is usually written 

^ = r,^ (16) 
p p 

where Fi is the adiabatic index (which need not be related to the various back- 
ground quantities). In terms of the Eulerian perturbations, this corresponds 
to 

5p Sp 



— = ^l^ + ^l^ 
p p 



Vlogp- -^Vlogp 
i 1 



(17) 
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At this point it is customary to introduce the so-called Schwarzschild discrimi- 
nant, As- It corresponds to the magnitude of the vector in the square bracket 
above, and since p and p depend only on a in our rotating model we have 

p da I ip da 

As we will see later, the Schwarschild discriminant plays a crucial role in stel- 
lar pulsation theory. Physically, radial variation in As corresponds to internal 
composition or temperature gradients in the star. In particular, the special case 
As = 0, which leads to the perturbations obeying the same equation of state as 
the background fluid (Fi = F), is often used for neutron stars. In the follow- 
ing, we will refer to such perturbations as being "isentropic" . The motivation 
for this assumption is that the temperature of most neutron stars is far below 
the relevant Fermi temperature, and they can consequently be considered as 
essentially having zero temperature! 

2.2. Rotational effects on pulsation modes 

We are now well equipped to discuss the effect that rotation has on the 
various pulsation modes of a star. To understand this issue it is helpful to 
recall how a mode calculation proceeds in the non-rotating case. For spherically 
symmetric stars one can separate the pulsations into two general classes. We will 
refer to these as spheroidal and toroidal perturbations. In relativistic studies 
these two classes are often called polar and axial perturbations (or even and 
odd parity perturbations), respectively. 

Spheroidal modes have displacement vectors of form 



e V-/o rr . Hi 



Sim, Himdg, ^^d^ ] Yim (19) 

ha 

where Yim{0,'^) are the standard spherical harmonics. Spheroidal perturba- 
tions are accompanied by variations in the pressure and the density which can, 
since they are scalar quantities, always be expanded in terms of the spherical 
harmonics. Hence, we have 

Sp^Y. ^PlmYlm (20) 
Im 

and similar for Sp. By analyzing the equations for a spheroidal perturbation one 
can deduce that a simple stellar model allows two distinct classes of pulsation 
modes (cf. the monograph by Unno et al.l for details). These have frequencies 

*However, it is worth pointing out that reaUstic neutron stars may still have As 7^ becapsp 
of radial variation in the chemical composition (say, a varying proton fraction in the core)E3. 
In view of this, whenever we refer to the "isentropic" case in this article, we intend this to 
mean stars with no internal entropy or composition gradients. 
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.1(1 + l)cl 



that are governed by 

« ^^^^ (21) 
where — Ap/Ap is the sound speed, and 

w -gS,rAs (22) 

where g = \dp/dr\/p is the local gravitational acceleration. These classes of 
modes are known as the p- and the g-modes, respectively. 
Toroidal modes have eigenvectors that can be written 

trv sm0 ; tr V^(^+l) 

where are the magnetic multipoles introduced by Thorne0. A non-rotating 
perfect fluid model has toroidal modes with non-zero freq uency , but in a model 
with a solid crust there are distinct toroidal shear modesEZlES. It should also 
be mentkmed that the relativistic theory yields both toroidal and spheroidal 
ly-modescil. 

Because of the symmetry of the non-rotating problem, modes corresponding 
to different I and m decouple. In fact, it is sufficient to consider the m = case 
since the non-zero m solutions then follow after a trivial rotation. Hence, there 
is no need to sum over the various I and m in the case of perturbed spherical 
stars. The case of rotating stars is much more complicated. First of all, the 
symmetry is broken in such a way that the various —l<m<l modes become 
distinct. As a first approximation one finds thatQ 

LOiin) = LOrin = 0) - mQ{l + Clrn) + 0{n^) (24) 

according to an inertial observer. Here, Cim is a function that depends on the 
mode-eigenfunction in a non-rotating star. Secondly, rotation has a significant 
effect on the eigenfunctions by coupling the various Z-multipoles. As the rotation 
rate is increased an increasing number of Yim^s are needed to describe a mode. 
One must also account for coupling between the spheroidal and toroidal vectors. 
To illustrate this, let us consider a toroidal mode that corresponds to a multipole 
m] in the = limit. After including the first rotational correction the 
eigenfunction of this mode becomes 



— — I 0, . „ du,, —Timde] Yi„i + ( Si^m, H^mde, —. — -du,]Y^rn (25) 
a \ sm I ^-^ \ sm I 



In other words, the [I, m] toroidal mode couples to the 1, m] spheroidal ones. 
The situation is analogous for rotationally modified spheroidal modes. 



2.3. Low-frequency modes 
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A rotating star has tw o sets of low- frequency modes, the spheroidal g-modes 
and the toroidal r-rnodestSnjI^'EH. To understand the origin and nature 
of these modes it is helpful to digress on the perturbations of non-rotating 
(spherical) stars. A spherical star has a degenerate spectrum of zero-frequency 
modes. These modes correspond to neutral convective currents, and the nature 
of the corresponding fluid perturbations can be understood in the following way: 
As follows from (|l^ and ([l3| ) stationary, isentropic, perturbations of a spherical 
star arc solutions to 

V • (pdu) = (26) 

and 

^Vp - -VSp ~ V(5$ = (27) 
P P 

Since these two equations decouple we can consider any solution as a superpo- 
sition of two distinct classes of solutionst3: 

5u ^ , 5p^5p = 5^ = Q 
Sit = , Sp, dp, (5$ nonzero 

Since any static self-gravitating perfect fluid must be spherical the second type 
of solution simply identifies a neighbouring equilibrium model. Hence, all sta- 
tionary nonradial perturbations of a spherical star must be of the first kind. 

What can we deduce about the fluid velocity for these stationary modes? 
First of all, one can readily verify that all toroidal displacements satisfy (p6|). 
In other words, the function Tim is unconstrained. In the case of spheroidal 
perturbations the situation is slightly different, and one finds that the following 
relation between Sim and Him raust be satisfied; 

-^{pr^Sim)-l{l + l)pr^Him^O (28) 
ar 

Still, one of the two functions is left unspecified also in this case. 

Before moving on to rotating stars, we need to consider also the case of 
non-isentropic perturbations. Then ( p7| ) provides a relation between the radial 
component of ^ and Sp and Sp. This will obviously not affect the above conclu- 
sion for toroidal modes: There is still a space of zero-frequency toroidal modes 
in a non-isentropic star. But the stratification associated with a non-zero As 
breaks the degeneracy in the spheroidal case, and gives rise to distinct g-modes, 
cf. (|2|). 

Having understood the nature of the zero-frequency modes in the = 
limit, it is easy to see what happens in the rotating case. For spinning stars the 
degeneracy of the zero-frequency modes is broken. The Coriolis force provides 
a weak restoring force that gives the toroidal modes genuine dynamicsflfl This 
leads to the so-called r-modes. The main properties of the r-modes can be 
understood by carrying out a simple exercise. Assume that the mode is purely 
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toroidal to leading order also in the rotating case. Then the motion is essentially 
horizontal and and 5p are of higher order in f2. Now consider the radial 

component of the equation we get by taking the curl of (]T^); 

9t(V X 5u)r + 2V X (l1 X 5u) = dt{V x 6u)r + Suh ■ Vhi^nr) = (29) 

In this equation subscripts r and h represent radial and horizontal components, 
respectively, and we have neglected terms of order 17^. We now translate this 
result into the inertial frame and get 



dt 



(V xSu) + 2n =0 (30) 



Here we can identify the quantity in the square bracket as the radial component 
of the total vorticity. This result thus shows_that the radial component of the 
vorticity is conserved by the r-mode motioiu. Furthermore, by inserting the 
toroidal eigenvector in ( p9| ) we readily find a first approximation of the r-mode 
eigenfrequency 

'^^ " WTT) ■ 

At this level of approximation different radial shells in the star are effectively 
uncoupled. Hence, the radial dependency of the r-mode eigenfunctions remain 
undetermined. To determine the radial behaviour of the modes we need to 
extend the analysis to order fi^. The reason for this can be understood as 
follows: From the linearised equations we can see that our leading order toroidal 
perturbation will not couple to the density and pressure perturbations until 
at order fl^. Indeed, if we write down the equations that describe terms of 
order Q we find that they decouple completely from the leading order toroidal 
perturbation. To linear order in fl the r-mode assumption is therefore consistent 
with Sp = Sp = = £,r = 0. As far as the higher order corrections to 
the r-mode frequency are concerned, one can argue that the next non-trivial 
contribution must be of order fl'^. Essentially, the frequency correction should 
have the same effect relative to the leading order frequency irrespective of the 
sense of rotation of the star. It cannot be that the r-mode frequency is different 
in magnitude for positive and negative f2. If that were the case the frequency 
would depend on the orientation of the observer, and change if he/she (say) 
decided to stand on his/her head, which clearly cannot make sense. Thus we 
anticipate that the r-mode frequency can be written in terms of a dimensionless 
eigenvalue uj2 as 



2mD. 

UJr 



1(1 + 1) 



1- LU2- 



M 



(32) 



We also expect that the r-mode eigenfunctions will take the form (g5|) with 
Tim = 0(1) and Si±im ^ Hi±im 0{fl'^), and that the pressure and density 
variations will both be of order fi^. 
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Table 1. Calculated r-mode frequency corrections uJ2 for constant density stars {n = 0) and n = 1 
polytropes. Results are shown both for the full calculation (nC) and the Cowling approximation 
(C). The results are for isentropic perturbations and the fundamental I = m modes. 

1 L02{n = 0) 0^2 = 1) 
nC C nC C 

2 0.765 0.913 0.398 0.453 

3 0.797 0.844 0.427 0.443 

4 0.730 0.749 0.399 0.405 



Once we have cranked through the algebra, the following picture emerges: 
For adiabatic non-isentropic perturbations, we have a Sturm-Liouville problemEj 
for the function Tim- This means that there will be an infinite set of r-modes for 
each combination of I and m. These modes are suitably labeled by the number 
of radial nodes in the respective eigenfunctions. A word of caution is in order, 
however. It is clear that we can reliably determine only a limited number of r- 
modes within the slow-rotation scheme: At some point the frequency correction 
will become large, and the assumption that it is a small perturbation will be 
violated. From (^2|) we can see that the scheme breaks down unless 

« ^ (33) 

Interestingly, the r-mode results for the isentropic case are rather different. 
When As — ori£^ can determine the eigenfunctions already at lowest order 
in the calculationca. Formally the calculation is taken to order il^, but when 
As — one can combine the r and components of ( p!^ ) to provide an equation 
relating the leading order contributions to (^g and This yields a rather simple 
equation for Tim'- 



sin - ^ (a^Tim) cos 0P„„ - (34) 
at) da 



Using standard relations for the Legendre functions one can show that this 
implies that we must have Tim = unless I ~ m- In other words, in the 
isentropic case r-modes can only exist for I — m, and in that case the above 
equation readily yields solutions of form 

Tu ~ a'-^ . (35) 

If we proceed to a full order mode calculatio we find the 

frequency corrections listed in Table ||. 

How can we understand the radical difference between the isentropic and 
the non-isentropic cases? A more detailed study of the r-mode eigenfunctions 
in the non-isentropic case shows that the fundamental mode (which has no 
nodes in the radial eigenfunctions) is well described by (|5|) also when As^ 0. 
Furthermore, the eigenfrequency of this mode remains virtually unchangecS as 
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As 0. Meanwhile, the frequency corrections for all the other r-modes change 
dramatically. This suggests the following explanation: A non-isentropic star has 
an infinite number of r-modes for each I and m ^ 0. These modes are distinct 
because of the stratification associated with the radial variation of As ■ In this 
sense these modes are close relatives to the spheroidal g-modes. The single 
r-mode that remains for ^ = m in an isentropic star is likely associated with 
the fact that the star is a sphere and not a cylindeiEj (which is relevant since 
the Coriolis operator has cylindrical symmetry). This "stratification" leads to 
the existence of a unique mode also in the isentropic case and explains why 
this mode is weakly dependent on variations in As- Only one mystery remains: 
What happens to all the other r-modes as the star becomes isentropic? We will 
return to this question in section 7.1. 

Detailed calculations show that the r-modes are such that Ap << 5p also to 
order . This means that the modes are well confined to the potential surfaces 
of the rotating configuration. In principle, this means that one would expect 
these modes to be associated with small variations in the gravitational potential. 
Given this, and the fact that the modes are more or less localized in the low 
density surface region of the star, the Cowling approximation (wherein is ne- 
glected) should be accurate for the r-modes. Calculations have verified thaLlhe 
modes are obtained to within a few percent in the Cowling approximationr^rH. 
This is illustrated in Table |^. 

What do the r-mode results imply for rapidly spinning neutron stars? A 
stable star can never spin faster than the rate at which matter is ejected from 
the equator. This corresponds to the Kepler limit, which is well approximated 

by 

Qk-^V^ (36) 

where p represents the average density of the corresponding nonrotating star. 
We can compare this estimate to an empirical formula, based on fully relativistic 
calculations for stellar models using realistic equations of state, which suggests 
thatiS'i 

w 0.78v/^ (37) 

Typically, the Kepler limit corresponds to a rotation period in the rang 
2 ms. It has been shown that O(ri^) slow-rotation stellar models are accurate 
to within a few percent compared to full numerical solutions. Consequently, one 
would expect a slow-rotation mode-calculation to be reasonably accurate even 
for the fastest spinning neutron stars. From (^3|) we immediately see that an 
r-mode solution ought to be reliable as long as 

|w2| « 5A/i.4i??o^-3 ■ (38) 
We have parameterized this formula in terms of canonical neutron star values: 



Mi.4 = M/1.4M( 



(39) 
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Rio = i?/10kni (40) 
P_3 = P/1 ms (41) 

Given this and the data in Table |^ we see that the fundamental r-mode as 
calculated in the slow-rotation approximation may be reliable also for the fastest 
spinning pulsars. 

2.4. Digression: The fluid motion 

Before we turn our attention to the interplay between r-modes, gravitational 
radiation and viscosity, it is relevant to discuss various features of the modes in 
some more detail. 

An important concept in the study of pulsation modes of rotating stars is 
the "pattern speed" of the mode. Given that every mode is proportional to 
^i{mip+u)t) j,g^j-^ that, holding a and 9 fixed, surfaces of constant phase are 
described by 

mip + ut — constant (42) 
After differentiation this leads to 

? = = - (43) 

which defines the pattern speed a of the mode. 

Having defined this concept, we make two observations. First of all we see 
that the pattern speed for a typical r-mode is 

(«) 



1(1 + 1) 

according to an observer rotating with the star. On the other hand, an inertial 
observer would find 

= " 1(1 + 1) ■ 

That is, although the modes appear retrograde in the rotating system an inertial 
observer would view them as prograde. 

Another useful result concerns the p-modes. Calculations for nonrotating 
stars show that the frequency of these modes increases (for a given overtone) 
with m slower than linearly. This means that as we increase m, the pattern 
speed of the p-modes will decrease, cf. (^. As a consequence, even though 
the high order p-modes have arbitrarily large frequencies, one can always find 
a mode with a very small pattern speed. The implications of this result will 
become clear in the following section. 

Let us now return to the r-modes and focus on the detailed fluid motion. We 
will concentrate on the I — m r-modes of isentropic stars. Introducing a suitable 
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dimensionless amplitude a, we can write the r-mnde velocity field (measured by 
an observer that is co-rotating with the star) aalj 

u«al7i?(|)'l^fe-'-* (46) 

We can readily illustrate this mode by drawing this vector field on the two- 
sphere, cf. Figure and Figure 1 of Saio^. This gives us an idea of the overall 
mode pattern. Every half period of oscillation the motion changes direction, 
and as we have already deduced, the pattern drifts backwards with respect to 
the sense of rotation (according to an observer on the star). 



3.2- 
JH 

2. 

2.6- 
2.4^ 
2.2 
2- 
1.8 
07.6 
1.4- 
1.2 
1 

0.8- 
0.6 
0.4 
0.2- 




o 



o 



o 



o 



Fig. 1. The large left frame shows the velocity field of the I = m = 2 r-mode at an instant in 
time. The narrow frame on the right is a qualitative picture of the actual motion of individual 
fluid elements (to leading order in the slow-rotation expansion). 

The velocity field in Figure |l| does not, however, give us a good picture of 
the motion of the individual fluid elements. To understand the nature of the 
r-modes better we consider fluid elements on the surface of the star. We can 



deduce from (46) that each fluid element moves according to 



(5r « , (47) 

59 — — oc asin6'cos[TO(^ -I- Wr*] , (48) 

5ip = — cx Q!sin6'cos6'sin[m(^ + LJr^] ■ (49) 
R 



From these equations we see that 



^ cxa^ (50) 



sin^ 6 SIT? 9 cos^ 9 
In other words, to leading order each fluid element moves on an ellipse the 
size of which is determined by the mode-amplitude a. As — > or tt the 
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ellipses approach circles and the radius shrinks to zero. We can also infer that a 
fluid element on the equator only moves up and down in the 9 direction. These 
features are illustrated in Figure 0. This is obviously only a first approximation, 
and one would expect higher order terms (both in f2 and a) to be important for 
rapidly spinning stars. We will return to briefly discuss the fluid motion and 
possible higher order secular effects in Sections 7.3-7.4. 

It is also worth trying to gain some understanding of the mode-amplitude a. 
Since we are working within a perturbative approach we naturally require that 
a be small, but what exactly does this mean? A useful measure is provided by 
the radial displacement. We know that the radial displacement is of order H,^, 
and we can approximate it using ( ^ ) together with equations (12)-(14) from 
Kokkotas and StergioulasS. In particular, if we consider a fluid element located 
on the surface of the star, we flnd that 



This estimates the "height" of the surface waves induced by the r-modes. Later 
on, when we discuss the possible effect that the unstable r-modes may have on 
the spin of a newly born neutron star, we will consider values of a of the order 
unity. Given that a neutron star has a radius of roughly 10 km, such values of 
a are obviously by no means "small" , as the correspond to waves of amplitude 
700 m on the surface of a star spinning at the Kepler limit. This is a very 
important insight, since it means that our discussion must be considered as (at 
best) an extrapolation into a regime where nonlinear effects should be relevant. 

Finally, it is worth emphasising that even though the r-modes are relatively 
exotic to astrophysicists they are very familiar to oceanographers and meteorol- 
ogists. The r-modes are analogous to the so-called Rossby waves in the Earth's 
atmosphere and oceans. These waves were first discussed by the Swedish mete- 
orologist Carl-Gustaf Rossby, and it is the association with Rossby waves that 
gives the r-modes their name. There are many interesting discussions of terres- 
trial Rossby waves in the literature. We should try to learn as much as possible 
from our colleagues in, for example, geop hys ics. Good starting points are the 
monographs by Greenspaneil and Pedloskyo. 

3. Gravitational-wave instabilities in rotating stars 

3.1. The CFS mechanism 

Many different instabilities can operate in a neutron star. The most familiar 
instability is due to the existence of a maximum mass (beyond which the star 
must collapse to form a black hole). Furthermore, above there are no 
equilibrium states (unless the star rotates differentially). Thus the star must 
shed some of its angular momentum before it can settle down into a stationary 




(51) 
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state. There have been many studies of this process in the hteraturalj'OE3. The 
so-called bar-mode instability grows essentially on the dynamical timescale of 
the star, and leads to gravitational waves that may well prove to be observable. 

It has long been recognized that the addition of various dissipation mecha- 
nisms may lead to secular instabilities, proceeding on a relatively long timescale 
depending on the magnitude of the added effect. The archetypal secular insta- 
bility is associated with viscosity. This instability was first studied by Roberts 
and Stewartson*^ (for more recent results, see Lindblom0). Essentially, vis- 
cosity admits a transition to a lower energy state by violating the conservation 
of circulation (which holds in a perfect fluid star). In a similar way, rotating 
neutron stars are generically unstable due to the emission of gravitational waves 
(which leads to angular momentum of the star not being conserved). That such 
an instability would operate was first established by Chandrasekhar@ for the 
Maclaurin spheroids. This intere£tiiLg_Lesult was subsequently put on a rigorous 



reiatmsjie 

tzBEae, 



footing by Friedman and SchutzcaEliO, who also proved that the instability is 
generic: All rotating perfect fluid neutron stars are unstable! 

The mechajiiain for gravitational-wave instability can be understood in the 
following wayeBES. Consider first a non- rotating star. Then the mode-problem 
leads to eigenvalues for w^, cf. Unno et alEJ, which in turn gives equal val- 
ues ±|a;| for the forwards and backwards propagating modes (corresponding to 
m = ±|m|). These two branches of modes are affected by rotation in different 
ways, cf. (p^. A backwards moving mode will be dragged forwards by the stel- 
lar rotation, and if the star spins sufficiently fast the mode will move forwards 
with respect to the inertial frame. Meanwhile, the mode is still moving back- 
wards in the rotating frame. The gravitational waves from such a mode carry 
positive angular momentum away from the star, but since the perturbed fluid 
actually rotates slower than it would in absence of the perturbation the angular 
momentum of the retrograde mode is negative. The emission of gravitational 
waves consequently makes the angular momentum of the mode increasingly 
negative and leads to an instability. This class of frame-dragging instabilities 
is usually refered to as Chandrasekhar-Friedman-Schutz (CFS) instabilities. It 
is easy to see that the CFS mechanism is not unique to gravitational radiation. 
Any radiative mechanism will do, and even though this possibility has not yet 
attracted much attention, one wo uld expect a similar instability to exist also 
for electromagnetic waves (see E£iJ and comments in section 7.3). 

The fact that the emission of gravitational radiation causes a growth in the 
mode energy in the rotating frame, despite the decrease in the inertial frame 
energy, may at first seem a bit strange. However, it can be understood from the 
relation between the two energies: 

Er ^ E,-nj . (52) 



From this we see that Er may increase if both Ei and J decreaseoEB. In other 
words, when the mode radiates away angular momentum the star can find a 
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Fig. 2. A schematic illustration of the conditions under which the CFS instability is operatin 
A perturbed star can be viewed as a superposition of a uniformly rotating background and 
nonaxisymmetric perturbation. A mode is unstable if it is retrograde according to an observer 
the fluid (left), but appears prograde in the inertial frame (right). 



rotational state of lower angular momentum and lower energy. Under these 
conditions the mode amplitude may grow. 

Let us take a closer look at some of the relevant ideas for this mechanism. 
The Friedman-Schutz criterion for instability relies on the so-called canonical 
energy Ec being negative. The canonical energy is defined as 

Be = \ j [pm\^ - p\u ■ ve? + tipi V • eT + f • vpv • c+ 

+ f.VpV-f +e*e*(V^V,p + pV,V,$)--^|V5$|2 (53) 

47rG J 

which is a conserved quantity. 

Similarly, we can define another conserved quantity 

Jc--Re j p^^C{^ti^+u■Vi^)dV (54) 

i.e. a canonical angular momentum. If Ec is negative at the outset and the 
system (the star) is coupled to another system (the radiation) in such a way Ec 
must decrease with time, then the absolute value of Ec will increase and the 
associated mode is unstable. Generally, an instability can be distinguished in a 
mode-independent way by constructing (canonical) initial data dt^] such that 
Ec is negative. However, to do this is not at all trivial. From the computational 
point of view the most straightforward task is to find the so-called neutral 
modes (which have zero frequency in the inertial frame), that signal the onset 
of instability. Furthermore, if we manage to proceed further and calculate a 
general mode of a rotating star, we can readily evaluate Ec to assess its stability. 
As was shown by Friedman and SchutzoE^, all mode solutions can be viewed 
as canonical data, and therefore it is sufficient to show that the displacement 
vector associated with the mode leads to i?c < to demonstrate the presence 
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of an instability. One can show that the condition Ec < is equivalent to the 
simple notion of a retrograde mode being dragged forwards, i.e. a change of sign 
in the pattern speed as viewed in the inertial frame. This important conclusion 
follows immediately from the relation 

Ec = - — Jc ^ (JiJc (55) 
m 

which is a general property of linear waves. Clearly, Ec changes sign when the 
pattern speed ai passes through zero. 

It is interesting at this point to contrast the gravitational- wave driven insta- 
bility to that due to viscosity. For uniformly rotating stars one can show that 
the combination 

SE - fl5J ^ Ec- nJc = - — Jc ^ (TrJc ^EcR (56) 

TO 

relating the first order changes in the kinetic energy and angular momentum to a 
mode-solution, is gauge- invariant. Ec^r can be viewed as the canonical energy in 
the rotating reference frame. Viscosity leads to Ec.r being a decreasing function 
of time. Comparing (^5|) to (|56| ) we see that the onset of the viscosity driven 
instability is signalled by the vanishing of the pattern speed in the rotating 
frame (dr = 0). 

3.2. Estimates for the f-mode instability 

Many families of neutron star oscillation modes contain members that, for 
a sufficiently high rotation rate, satisfy the CFS instability criterion. We have 
already pointed out that the large to p-modes may have arbitrarily small pattern 
speeds in a spherical star. Combining this result with the qualitative formula 
( p^ we can deduce that these p-modes will become unstable already at relatively 
slow rotation rates. Because it has no nodes in its radial eigenfunctions the 
fundamental p-mode (the /-mode) radiates gravitational waves more efficiently 
than the various overtones, and consequently leads to the strongest instability 
for any given to. 

So why is it that, if all rotating stars are generically unstable, we observe 
millisecond pulsars? This is a crucial question, the answer to which follows 
from several important results. First of all, we have seen that the large m 
modes become unstable at the lowest rates of rotation. But as Comina£3S have 
shown (for the Maclaurin spheroids) the growth rate of these modes decreases 
exponentially with to. In order to play an astrophysical role a mode must grow 
fast compared to the evolutionary timescale of the star. In practice, this means 
that one would not expect modes with m significantly larger than (say) 10 to 
be relevant. Secondly, it turns out that gravitational radiation and viscosity 
compete. Viscosity tends to suppress the CFS instability, and since it operates 
locally on short length scales it is more effective for high to modes. Detailed 
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calculationsOEJ have shown that viscosity suppresses the CFS instabihty in the 
p- modes for m > 5. 

To estabhsh a rehable consensus regarding the stabihty properties of rotating 
neutron stars is a high priority issue. However, despite a considerable amount 
of work being devoted to the task this subject may still hide many of its secrets. 
Given that a general criterion for stability is outstanding, we must attempt 
to calculate all the various modes of a rotating star to see if any of them are 
unstable. This is a formidable task, but considerable progress has been made. 
Currently, most of our understanding of unstable modes comes from Newtonian 
calculations, in particular for the Maclaurin spheroids. 

Until recently, the conventional wisdow was that the //p- modes would lead 
to the strongest gravitational-wave instability. This notion dates back_tp the 
very first investigations of the CFS mechanism. Friedman and SchutiSl sug- 
gested that "we can probably safely conjecture that thep-modes of nonisentropic 
stars are the most important ones for stability, since they typically involve larger 
density changes than the g-modes do, and so will radiate gravitational waves 
more effectively" . This would seem to be a reasonable assumption, and it ex- 
plains why studies of the CFS instability in neutron stars remained focussed on 
the //p-modes for almost twenty years. 

Before proceeding to discuss the new ideas regarding the CFS-instability 
that have emerged in the last couple of years, we will summarise the status of 
the field as of a few years ago. The current understanding of the instability 
in the /-modes of a rotating neutron star is essentially based on two different 
bodies of work. 

Newtonian results for the /-modes of rapidly rotating stars can be used 
to estimate the competing influences of gravitational radiation reaction and 
various viscosities. This leads to a prediction that the m ~ A mode provides 
the most stringent limit on rotation, and that that the CFS instability would 
limit the rotation of a normal fluid star toE3 > O.OSrii^. Similar estimates 
for superfluid stars, suggest that dissipation due to the so-called mutual friction 
(see section 4.5) will completely suppress the instabilityEj. 

A fully relativistic calculation of pulsation modes of rotating stars is still 
outstanding. The best results so far concern neutral (zero frenuency in the 
inertial frame) modes that would signal the onset of instabilityo, and modes 
calculated within the Cowling approximationElH. As could perhaps have been 
anticipated, these results indicate that relativistic effects strengthen the insta- 
bility. For example, while Newtonian calculations predicted that the m = 2 
/-mode does not become unstable for $7 < Q.k the relativistic results show that 
this mode can in fact become unstable and that it may provide the strongest 
constraint on the rotation rate of a relativistic star. 

4. The r-mode instability 



Andersson & Kokkotas 19 



The discovery that the r-modes of a rotating neutron star are generically 
unstable due to the emission of gravitational wavesQQ may have come as a 
slight surprise, but in retrospect the result is rather obvious. The real surprise 
is that this was not realized earlier. After all, the r-modes have been discussed 
in Newtonian stellar pulsation theory for the last tsrenty years or so (following 
a pioneering 1978 paper by Papaloizou and PringleEa), and their properties are 
described in monographs on the subjectQ. Since the modes are prograde in the 
inertial frame and retrograde in the co-rotating frame, cf. (^4|)-(^5|), they satisfy 
the Newtonian criterion for the CFS-instability. Hence, one can easily deduce 
that the r-modes ought to be unstable. That this instabilitv is generic also in 
the relativistic case was first shown by Friedman and Morsinm, who proved that 
there exist toroidal initial data for which the canonical energy Ec is negative. 

In perfect fluid stars the r-modes are unstable at all rates of rotationQ. This 
is an important observation, because it indicates that the instability may well be 
astrophysically relevant even though the r-modes do not lead to large variations 
in the density and therefore would not intuitively be associated with strong 
gravitational waves. The reason one might expect the r-mode instability to be 
competitive with the more familiar instability in the /-modes is that an /-mode 
has nonzero frequency in a spherical star. Thus, the star must be spun up to 
a critical rotation rate (flc) in order for the /-mode to go unstable, cf. (p^. 
Because the r-modes are unstable at all rates of rotation 

their growth times may become rather short already at spin rates below fie, 
which could then lead to the r-mode instability being stronger tJian the /-mode 
one. The first attempt to calculate the r-mode growth timeda suggested that 
this would be the case. The growth times for the r-modes were estimated to be 
very short indeed. It was suggested that the amplitude of an unstable r-mode 
would grow by a factor of two in a second for a typical neutron star spinning with 
a period of 1.6 ms. However, these estimates came with a serious disclaimer: 
The method used to determine the modes, and infer the growth times, was not 
designed to study low-frequency modes. 

Nevertheless, these first estimates provided a strong motivation for more 
detailed work. As a first step towards understanding this new instability better 
it seemed sensible to obtain estimates that allowed an immediate comparison 
with the established /-mode results (recall that these estimates suggest that 
the /-mode instability becomes active at rotation rates above 0.95ri/f). Such 
a comparison would establish whether the r-mode instability is of relevance 
astrophysically or if it is just a peculiarity associated with general relativity. 
To make the desired comparison one must calculate the characteristic growth 
time of the modes, as well as the damping times due to various dissipation 
mechanisms that may be relevant. As in the case of the /-mode, the stability 
of the r-modes depends on the competing influences of gravitational radiation 
reaction (that drives the mode) and viscous damping. In order to establish the 
astrophysical relevance of the instability we must confirm two things: First, the 
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unstable mode must grow on a timescale that is astrophysically "short", for 
example, significantly below the age of the universe. Secondly, the instability 
must be able to win the tug-of-war against all relevant dissipation mechanisms. 

In this section we will discuss various estimates that have been made, espe- 
cially for damping mechanisms that may counteract the growth of an unstable 
r-mode. The estimates we present are made in the framework of Newtonian 
gravity. Since the r-mode instability is a truly relativistic effect this obviously 
makes these results somewhat ad hoc. But we must remember that there are 
as yet no fully relativistic calculations of pulsation modes of rapidly rotating 
neutron stars. In absence of such results we expect Newtonian estimates to pro- 
vide useful indications. One would not expect more detailed models, eg. in full 
general relativity, to completely change the Newtonian picture (we comment on 
some recent results regarding relativistic r-modes in section 7.2). After all, the 
relativistic corrections are likely to be of the order of (say) 20% and there are 
far larger uncertainties associated with other aspects of this problem (eg. the 
microphysics). 

We will adopt the standard strategy for estimating instability timescaleB. 
First we assume that the true mode-solution is well represented by the solu- 
tion to the non-dissipative perturbation equations (as described in Section 2.3). 
Then we use these solutions to evaluate the effect of the various dissipation 
mechanisms and add their respective contributions to the rate of change of the 
mode energy dE/dt. Finally, we can verify that the first assumption is justi- 
fied by checking that the estimated growth/damping time is considerably longer 
than the oscillation period of the mode. 

We estimate the timescale td of each dissipation mechanism by assuming 
that the eigenfunctions are proportional to exp(i/t(j). Then we recall that the 
mode-energy follows from the square of the perturbation and deduce that 



dE _ 2E 
dt td 

where 



(57) 



E^IJ pM'dV « ^l±^u;'^ j\a'\Tu\'da (58) 



is the energy of an r-mode measured in the rotating frame. Using (|4^) we can 
deduce that (for an n 1 polytrope) 

E w l0^^a^Mi,iRloPl3 erg (59) 

We focus our discussion on the I — m modes because they have been shown 
to lead to the strongest instabilityil'0. Note also that we only keep the leading 
order terms in (|58| ) and the dissipation integrals below. For simplicity, we 
also use the Cowling approximation, an assumption that only introduces an 
inaccuracy of a few percentSS, cf. Table 0. This is rather good precision in 
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an astrophysical context, and there are far larger uncertainties in our current 
understanding of the r-mode instabihty. 

IdeaUy one would like to address issues regarding the relevance of the supranu- 
clear equation of state for an active mode-instability. However, to use a realistic 
equation of state in a Newtonian calculation is not particularly meaningful. The 
reason for this is that there is no one-to-one correspondence between Newtonian 
and relativistic stellar models in the sense that, for a given equation of state, 
the same central density leads to stars with comparable masses and radii in the 
two theories. Thus we believe that a calculation for realistic equations of state 
in the Newtonian framework can be misleading. Still, it is interesting to try 
to understand whether the overall properties of the equation of state, eg. the 
stiffness, affect the instability. This issue can be studied by considering differ- 
ent polytropic models, and in the following we consider results both for n = 1 
polytropes and constant density stars (represented by the limit n — > 0). These 
results ought to be indicative for realistic neutron stars given that the average 
adiabatic index of most realistic equations of state lies in the range n ~ — I. 

4.1. Gravitational radiation reaction 

The first step in establishing the relevance of the r-niode instability cor- 
responds to showing that the modes grow on an astrophysically interesting 
timescale. To estimate the gravitational-radiation reaction we use the standard 
post-Newtonian multipole-formulas (see ThorneE3). The gravitational- wave lu- 
minosity associated with a pulsation mode (measured in the rotating frame) is 
then estimated using 



dE 
'dt 



^-u;rJ2^i^?^H\^Di„,f + \SJi,nf) , (60) 
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where 



AT - V ^ -^A-- ^7 (f..^ 



AttG {l + l){l + 2) 
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The first term in the bracket of ( |60| ) represents radiation due to the mass mul 
tipoles. These are determined by 



SDi,„. = J Spa%*„^dV . (62) 

The second term in the bracket of (|60| ) corresponds to the current multipoles, 
which follow from 



SJhn = -\ I a'ipSu + dpn)Yi'^*dV , (63) 



where are the magnetic multipole: 
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Table 2. Estimated growth times due to current and mass-multipole gravitational radiation. The 
results are for the single I = m r-mode that exists for isentropic stars, and are given both for 
constant density stars (n = 0) and n = 1 polytropes. We do not show mass multipole results for 
the p = constant case, since (5p = in that case. 
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An order count based on the above formulas suggests an interesting result: 
For I = TO r-modes, the dominant contribution to the gravitational radiation 
comes from the first term in (^3|). That this is the case can be seen as follows. 
Recall that the I = m modes have toroidal displacement to leading order, while 
the spheroidal perturbations enter at order fl^. Furthermore, we know that if the 
toroidal component corresponds to the /th multipole, the spheroidal components 
will correspond to / -f 1. This means that, as a gravitational- wave source the r- 
modes are quite unusual. Since they are primarily perturbations of the velocity 
field in the star, with little disturbance in the star's density, the gravitational 
radiation that they emit comes primarily from the time-dependent mass currents 
(we see that SDim ^ while 6 Jim ^ This is the gravitational analogue of 
magnetic multipole radiation. In fact, the r-mode instability is unique among 
expected astrophysical sources of gravitational radiation in radiating primarily 
by gravitomagnetic effects. 

After inserting the leading order r-mode eigenfunction in the relevant current 
multipole term, we find 
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(64) 



for the rate of energy loss due to gravitational waves. This then leads to an 
estimated growth timescale for the instability 



-tr 



7i,r— 1 p — 2/ p2/ + 2 



(65) 



where the negative sign indicates that the mode is unstable. 

The estimated current multipole timescales are given (for I = m r-modes of 
an isentropic star) in Table B. Th e l isted results are consistent with original ones 
obtained by several authorslMoEJCj From the tabulated results we can deduce 
two things. First of all we see that the timescale increases by roughly one order 
of magnitude with each Thus, higher multipoles lead to significantly weaker 
instabilities and the / = to = 2 r-mode will be the most important. Secondly, we 
can see that the results for constant density stars and n = 1 polytropes differ by 
roughly a factor of two. This provides a useful illustration of the uncertainties 
associated with the supranuclear equation of state. 
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It is worth pointing out that the conclusion that the current multipoles 
dominate the gravitational-wave emission is not necessarily true for all the r- 
modes. It is certainly true for isentropic stars, since such stars have no "pure" 
r-modes for I 7^ m. But in the general non-isentropic case we have to consider 
the infinite set of modes that exist for all m ^ 0. An order count in ( |60[ ) 
immediately shows that the situation is not so simple for the I ^ m modes. In 
general, the mode-solution for these modes have contributions to the toroidal 
eigenfunction described by the [l,m] multipole, but the spheroidal part and the 
density perturbation (proportional to fl^) correspond to [/±1,to]. Consequently, 
the current multipoles and the mass multipoles contribute to ( |60|) at the same 
order in for many of the non-isentropic modes. Thus mass multipole radiation 
could play an important role. It is also relevant to stress that the mass multipoles 
are far from irrelevant for the isentropic I = m modes. From ( |62|) wc find that 
the leading contribution to the mass multipole radiation follows from 
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where we have introduced 



This leads to 



5pim = pgaQn 



(67) 



where the estimated values for tmm are listed in Table Here it is interesting 
to note that the mass and current multipole results scale with the size of the 
star in rather different ways. This means that the relative importance of the 
two is significantly different for stars of different size. 



4.2. Viscous damping 

The estimates given in Table ^ show that the r-modes grow rapidly enough 
to be of potential significance. Still, they must also overcome damping effects 
due to different physical mechanisms, most of which are poorly understood and 
therefore difficult to model in a realistic fashion. 

As in the case of the unstable /-modes the main dissipation in a hot newly 
born neutron star is due to viscosity. For the r-modes to be relevant they must 
grow fast enough that they are not completely damped out by viscosity. To 
assess the strength of the viscous damping of the r-modes we implement the ap- 
proximations that have been used to study unstable /-modes in the Newtonian 
context0. For the simplest neutron star models two kinds of viscosity, bulk 
and shear viscosity, are normally considered. These are due to rather different 
physical mechanisms and we discuss them separately below. 



4.2.1. Shear viscosity 
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At relatively low temperatures (below a few times 10^ K) the main vis- 
cous dissipation mechanism in a fluid star arises from momentum transport 
due to particle scattering. In the standard approach these scattering events 
are modelled in terms of a macrosopic shear viscosity. In a normal fluid star 
neutron-neutron scattering provides the most important contribution. In a su- 
perfluid the situation is a bit more subtle. Electron-electron scattering leads to 
the dominant contribution to the shear viscosity, but one must also account for 
exotic effects like the scattering off of vortices in the superfluid, cf. Section 4.5. 

The effect of shear viscosity on the r-modes can be estimated fromE 



dE 
'dt 

where the shear ij„h follows from 



= -2j rjSa^'Sal.dV (69) 
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After working out the angular integrals we are left with 
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(71) 

Above the transition temperature at which the neutron star becomes super- 
fluid (several times 10^ K), the appropriate viscosity coefhcient (due to neutron- 
neutron scattering) is 

r;„„ = 2 X lO^V^s^rg-^g/cms . (72) 

where 

Pi5 - p/lQiSg/cm^ (73) 
Tg = r/lO^K (74) 

This leads to an estimated dissipation time-scale due to shear viscositySSS 

t,v~A?i.4 ^10 ^9 j 6.7x lO^s n = l ^^^> 

To obtain these results the star has been assumed to be isothermal. This should 
be a reasonable approximation apart from for the first few moments of a neutron 
stars life. These estimates of the shear viscosity are expected to be relevant for 
the first months of the life of a hot young neutron star. 

Once the star has cooled below the superfluid transition temperature, the 
above viscosity coefficient must be replaced by 

ryee = 6 X 10iV?5?V'g/cnis . (76) 
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which follows from an analysis of electron-electron scattering. This leads tol 

isv ~ MiA^ioh \ 2.2 X 10^ s n = 1 

The above shear viscosity estimates were based on standard Navier-Stokes 
theory as developed by, for example, Landau and LifshitzEll. However, this 
theory is somewhat pathological in that it allows dissipative signals to travel 
faster than light. A causal theoryof dissipation has been developed following 
initial work by Israel and StewartEj. This alternative description differs from 
Navier-Stokes theory in some potentially important ways. The main feature of 
the Israel-Stewart description is that there is a coupling between rotation and 
the gradients of both temperature and momentum. Intuitively one might expect 
these additional effects to be small, but there are two situations where they may 
be relevant. For high frequency oscillations the large time-derivatives may lead 
to significant corrections to the Navier-Stokes results. Similarly, rapid rotation 
can lead to a strong coupling. The latter effect could certainly be relevant for 
the r-modes (or indeed any other mode Jn a rapidly spinning star). Inspired 
by this possibility, Rezania and Maartena^ have estimated the relevance of the 
coupling between vorticity and shear viscosity for the r-modes. Their results 
show that the correction due to the vorticity coupling can be large, especially 
for low temperatures. 

4.2.2. Bulk viscosity 

At high temperature (above a few times 10^ K) bulk viscosity is the dominant 
dissipation mechanism. Bulk viscosity arises because the pressure and density 
variations associated with the mode oscillation drive the fluid away from beta 
equilibrium. It corresponds to an estimate of the extent to which energy is 
dissipated from the fluid motion as the weak interaction tries to re-establish 
equilibrium. The mode energy lost through bulk viscosity is carried away by 
neutrinos. 

In the previous section we found that we could calculate the shear viscosity 
timescale from the leading order contribution to the r-mode fluid motion, cf. 
( [n]) . This is not the case for the bulk viscosity. To assess its relevance we need 
the Lagrangian density/pressure perturbation. But for r-modes these perturba- 
tions vanish to leading order in rotation. This means that the mode calculation 
must be carried at least to order fl^ if we want to estimate the bulk viscosity 
timescale. In particular, we need 

dE 



dt 



= -J(:M' (78) 



bv 

where Sa is the expansion associated with the mode, defined by 
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and 

Ap = Sp+ = Pga J2 [Om - Sim] Yira , (80) 

Im 

where we have used the fact that dp/da — —pg for a Newtonian modeL The 
detailed calculation of the relevant eigenfunctions has been described by, for 
example, Saio0. 

After doing the angular integrals we get (for I = m) 



dE 
lit 



bv 



C(^] [Q+ii-Si+ufda. (81) 



In the case when /3-equilibrium is regulated by the so-called modified URCA 
reactions the relevant viscosity coefficient is 

C^j, = 6 X lQ^'>p\,Tl (^) g/cms . (82) 

This leads to an estimated bulk viscosity timescale for n = 1 polytropes0 

tbv ~ 2.7 X 10"Mi.4i?ro^P^3rg-'^ s (83) 

This result, which was obtained within the Cowling approximatioii agrees well 
with a calculation including the perturbed gravitational potentialcj. 

We note that it has been speculated that the proton fraction may be large 
enough to make direct URCA reactions possible in the cores of neutron stars. 
If this is the case, the bulk viscosity coefficient becomes significantly larger 
than (|8^). However, direct URCAnrocesses will only be relevant if the proton 
fraction is above several percentElo. Since this is unlikely to be the case in the 
surface regions where the r-modes are mainly located we will not consider the 
direct URCA bulk viscosity in detail here. 

Let us also point out that, since very hot stars are no longer transparent to 
neutrinos, the above estirnates for the bulk viscosity will not remain relevant at 
very high temperatures lj. This may be of some relevance for the very early 
phase of r-mode growth in a newly born neutron star. 

Another point, tha.t may turn out to be of crucial importance, was recently 
raised by Jones EZI. He pointed out that the bulk viscosity result would be 
significantly different if the presence of hyperons in the neutron star core was 
accounted for. Not only is the associated viscosity coefficient stronger than as- 
sumed above, the temperature dependence is also different in the hyperon case 
(the coefficient scales as rather than T^). This makes the hyperon bulk 
viscosity relevant at low temperatures. Jones argues that the r-mode instability 
is almost completely suppressed in the case when hyperons are present through- 
out the entire star. This is, however, not a particularly realistic assumption. 
In order to assess the extent to which hyperons affect the estimated instability 
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time-scales calculations that allow for the presence of exotic particles (as pre- 
dicted by modern equation of state results) need to be performed. At the time 
of writing, no such results are available. 

Since there has been some debate in the literature regarding the estimated 
timescales for the bulk viscosity it is relevant that we try to identify the reasons 
why various given estimates differ. Fipt of all, let us discuss the constant density 
estimate of Kokkotas and StergioulasEj. As is easy to see from (^9|), one cannot 
rigorously calculate the expansion associated with a mode in a constant density 
star. After all, the assumption p — constant immediately leads to Ap ~ 0. Still, 
one can estimate of the relevant effect (following Cutler and Lindblom0) in the 
following way: Assume that the Lagrangian variation in the pressure takes the 
following form 

^.f(^)'.-«„ m 

where e is a dimensionless normalisation constant. This formula then leads to 
the anticipated behaviour as F — > oo, i.e. as p — > constant, and the formula also 
agrees with the result for /-mode oscillations^ in compressible stars (to within 
a factor of two). Kokkotas and Stergioulasc2l use this approximation also for 
the r-modes, and arrive at an estimated bulk viscosity timescale (actually using 
r = 5 in (H)) 

tbv ~ 2.4 X 10^° M^-jRloTg^ Pis s (85) 

Compared to ( p3| ) this is an underestimate of the strength of bulk viscosity 
dissipation by roughly one order of magnitude. 

An estimate that was equally far away from the correct value was given by 
Lindblom, Owen and MorsinkE3. They tried to approximate the bulk viscosity 
using only first order (in Q) results. In particular, they assumed that Ap w 6p. 
This assumption is, however, not warranted for the r-modes. As we discussed in 
Section 2.3, one of the fundamental r-mode properties is that the fluid essentially 
moves on isobars in the perturbed configuration. Hence, we will have \Ap\ << 
\Sp\, which results in Lindblom, Owen and Morsink overestimating the strength 
of the bulk viscosity by roughly one order of magnitude. 

Finally, the first published second order (in fl) results (due to Andersson, 
Kokkotas and Schuta^) also differ considerably from (|3|), essentially in the 
scaling with the rotation rate. This was due to a typographical error in the 
numerical code used to calculate the mode eigenf unctions. Once this mistake is 
corrected, we arrive at (^3|). 

4.3. The r-mode instability window 

Armed with the estimates obtained in the previous two sections we can ad- 
dress the main issue of interest: Will the r-mode instability be relevant for 
astrophysical neutron stars? First of all, it is easy to see that the r-modes will 
only be unstable in a certain range of temperatures. To have an instability we 
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need tgw to be smaller in magnitude than both tgv and tbv From the various 
estimates we immediately see that the dissipation due to shear viscosity kills 
the mode at low temperatures, while the bulk viscosity dominates at high tem- 
peratures. In fact, by comparing (^7|) to ( p5[ ) we deduce that shear viscosity will 
completely suppress the r-mode instability at core temperatures below 10^ K. 
Similarly, bulk viscosity will prevent the mode from growing in a star that is 
hotter than a few times 10^" K (but see the proviso in section 4.2.2 regarding 
the relevance of the bulk viscosity at high temperatures). In the intermediate 
range there is a window of opportunity where the growth time due to gravita- 
tional radiation is short enough to overcome the viscous damping and drive the 
r-mode unstable. 




Fig. 3. The critical rotation rates at which shear viscosity (at low temperatures) and bulk viscosity 
(at high temperatures) balance gravitational radiation reaction due to the r-mode current multi- 
pole. This leads to the notion of a "window" in which the r-mode instability is active. The data 
in the figure is for the I = m = 2 r-mode of a canonical neutron star [R = 10 km and M = IAMq 
and Kepler period P[f ^ 0.8 ms). 



This instability window is usually illustrated by the critical rotation period, 
above which the mode is unstable, as a function of temperatureEiS'El. We find 
the relevant critical rotation rate by solving for the zeros of 



2E dt tcrvi id: 



ISS 



for a range of temperatures. A typical result (for canonical neutron star param- 
eters) is shown in Figure ^ For canonical values, the Kepler limit corresponds 
to Pk — 0.8 ms. In the figure we show the critical period curve for an n = 1 
polytrope. From this data we can deduce that a rotating neutron star with a 
core temperature of 10^ K (i.e. a few months old) is unstable at rotation periods 
shorter than 25 ms. There is, of course, a large uncertainty associated with this 
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prediction. If wc, for example, assume that the instabihty sets in at roughly the 
same Pk/P for different equations of state (which is supported by the n = 
and 1 results), then the uncertainty in the Kepler period {Pk ~ 0.5 — 2 ms) 
for realistic equations of state suggests that the instability will be relevant at 
rotation rates faster than P « 10 — 40 ms at some core temperature. 

Figure || shows the different strengths of the normal fluid shear viscosity 
(due to neutron-neutron scattering) and the superfluid counterpart (due to 
electron-electron scattering). It is interesting to note that the superfluid is, 
rather counterintuitively, more dissipative that the normal fluid. We also indi- 
cate the strength of the / = m = 2 r-mode mass multipole. As a comparison 
we should recall that the m = 4 /-mode (believed to be the "most unstable" 
spheroidal mode) becomes unstable above Pk/P ~ 0.95. In other words, the 
r-mode instability provides a much stronger constraint on neutron star rota- 
tion than the instability in the /-mode does. In fact, even the comparatively 
weak radiation reaction associated with the r-mode mass multipoles leads to an 
instability that completely dominates the /-mode one. 




Fig. 4. The r-mode instability window for a canonical neutron star (the data is for an n = 1 
polytrope). The lower two curves illustrate the difference between shear viscosity due to neutron- 
neutron scattering (solid line) and electron-electron scattering (dashed curve). The latter should 
be relevant for a superfluid star (i.e. below a few times 10^ K). It should be noted that the 
superfluid shear viscosity is stronger than that of a normal fluid. We also show the instability 
window that arises if we only include the mass multipole radiation from the I = m = 2 r-mode. It 
is useful to compare these results to those for the m = 4 /-mode, that suggest that this mode is 
unstable for Pk/P > 0.95 or so. 

Finally, we illustrate the relative importance of the higher I — m r-modes 
in Figure ^ It should be noticed that, even though the growth time increases 
with roughly one order of magnitude for each / the higher multipole r-modes 
still provide a severe constraint on neutron star rotation rates. 

The results in Figures are encouraging and perhaps indicative of the 
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Fig. 5. The r-mode instability window for a canonical neutron star (the data is for an n = 1 
polytrope) and different Z-multipoles. 

astrophysical role of the r-mode instability! But more than anything else they 
emphasize the urgent need for an improved understanding of the various ele- 
ments in our model. We will now discuss various attempts to improve on the 
detailed physics included in the r-mode calculations. 

4.4. Superfluid mutual friction 

As a neutron star cools below a few times 10^ K the^xtreme density in the 
core is expected to lead to the formation of a superfluidOH. This means that, 
while it may not be wholly inappropriate to model a newly born neutron star 
as a simple "ball" of a more or less perfect fluid, a model of an old neutron 
star must certainly be more complicated. In many ways it is useful to think 
of an old (which in this context means: older than a few hours/days) neutron 
star as a "layer cake" . A simple, yet reasonable, picture of an old neutron star 
consists of i) a core in which superfluid neutrons (in a condensate) coexist 
with superconducting protons (above p w 1.5 x 10^'* g/cm"^), ii) a crust: a 
lattice of nuclei permeated by superfluid neutrons in a ^ 5*0 state (the superfluid 
reaches out to the neutron drip density ~ 10^^ g/cm^), and iii) a fluid ocean. 
A canonical neutron star would have a core of 9 km, a crust of 1 km and a thin 
ocean 50-100 meters deep. The superfluid constituents play a crucial role in 
determining the dynamical properties of a rotating neutron star. In particular, 
the interplay between the lattice nuclei and the superfluid in the inner crust 
is a key agent in the standard model for neutron star glitches. Intuitively, the 
presence of a superfluid should have a considerable effect also on the various 

tWe n ptc ii il l rFiassing that the r-modes may also be unstable in the fastest rotating white 
dwarfsllEj'Ei'Ea; the so-called DQ Herculis stars. 
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modes of pulsation. To some extent this is already indicated by the results in 
Figure ^, and when the presence of new dissipation mechanisms is accounted 
for it turns out that a superfluid may be more dissipative than a normal fluid! 
The most important new effect is the so-called mutual friction that arises from 
scattering of electrons off of the neutron vortices (recall that a superfluid mimics 
large scale rotation by forming a large number of vortices). This scattering is 
greatly enhanced by the strong magnetic field (due to the entrained protons) 
within the vortices. 

The motion associated with a neutron star pulsation mode is of large scale 
compared to both the neutron vortex and superconducting proton magnetic flux 
tube separations (~ fO^'^ cm and ^ 10^^" cm, respectively). This means that it 
is sufficent to include the averaged dynamics of many vortices and flux tubes in 
a model for superfluid stellar pulsation. This leads to models consisting of two 
distinct, but dynamically coupled, fluids representing the superfluid neutrons 
and the "protons" (a fluid that contains all charged components in the star). 
The dynamics can then be described in terms of two velocity fields, described 
by two sets of equations similar to ( [l^ ) and ( |l3|) with additional coupling terms 
(see equations (1) to (4) in Lindblom and MendellEl) . It is interesting to note 
that this siipprfliiid ynnHpl leads to a new set of pulsation modes already in the 
nonrotating caseOlIj. In order to discuss the difference between these modes and 
the familiar //p- modes it is useful to introduce a variable SP, which represents 
the "departure from /3-equilibrium" in the oscillating fluid. Then the standard 
fluid modes all have small S/3. For these modes the superfluid neutrons and the 
protons move more or less together. In contrast, the superfluid modes have (5/3 
large, and the neutrons and protons are largely countermoving. 

Superfluid mutual friction has been shown to cninpletely suppress the insta- 
bility associated with the /-mode in a rotating starE3. Plausibly, it is a dominant 
damping agent also on the unstable r-modes. We cariassess the relevance of mu- 
tual friction in a hand- waving way following MendellEU: For a longitudinal wave 
travelling perpendicular to the rotation axis the influence of mutual friction can 
be estimated as 

Here A is the wavelength of the oscillation (which we will take to be A ~ i? w 
10 km), and 

r„,f < 10i°T-2 cmVs (88) 
Comparing this to our estimated shear viscosity timescale we see that 

t^t > 200 s (89) 

and that this timescale is independent of the temperature. This is in clear 
contrast with the shear and bulk viscosities, both of which are sensitive to 
changes in the core temperature. This result suggests that, even though mutual 
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friction may be important it will not completely suppress the r-mode instability. 
In fact, we would expect a detailed calculation to lead to timescales considerably 
longer than this upper limit. This is since mutual friction would have the 
largest effect on waves that travel perpendicular to the rotation axis (as in 
our example) while the velocity field of a pulsation mode has a considerable 
component parallell to the axis. 

In the "macroscopic" picture, the strength of mutual friction can be ex- 
pressed in terms of the "entrainment parameter" 

e=^(!I^-l] (90) 

where pn and pp are the neutron and proton densities, respectively, and nip and 
m* are the bare and effective proton masses. Estimates of the ratio of these two 



masses suggest thatL3 



to: 



0.3<^<0.8 (91) 

TOp 

Steps towards more detailed studies of superfluid r-modes have been taken by 
Lindblom and MendellEa. They have calculated the modes in a superfluid star 
and shown that they are identical to those of a normal fluid star to lowest order 
in fl. Therefore the growth timescale due to gravitational radiation reaction 
remains essentially unchanged in a superfluid star, cf. (|6^). For the specific 
two-fluid stellar model considered by Lindblom and Mendell, the shear viscosity 
timescale turns out to be roughly a factor of 3 shorter than in the normal fluid 



case, cf. (75). As far as the mutual friction is concerned, the calculations of 
Lindblom and Mendell do not provide a complete answer to the question of 
whether this mechanism will suppress the r-mode instability or not. 

Lindblom and Mendell find that the typical result for r-mode dissipation 



due to mutual friction ii 



W ~ 2 X lO^Plg s . (92) 



This means that the r-niode instability window in (say) Figure |^ would be es- 
sentially unaffected by the inclusion of mutual friction. However, as manifested 
in Figure 6 of Lindblom and Mendell, the result is sensitive to changes in the 
entrainment parameter. In particular, there are critical values of e for which 
the mutual friction timescale becomes very short (and for which the unstable 
r-mode would be completely suppressed). It is interesting to note that 6(3 is a 
key parameter in deciding the importance of mutual friction. The critical val- 
ues of e for which the mutual friction timescales are short correspond to modes 
for which 5f3 is comparatively large. As already mentioned, the protons and 
neutrons are essentially countermoving for large Sp. Given that mutual friction 
tends to damp out the relative motion between the neutrons and the protons, 
it is natural that modes with large 6(3 should be strongly affected by it. The 
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results of Lindblom and Mendell seem to suggest that there are values of the 
entrainment parameter for which the r-modes are in "resonance" with a super- 
fluid mode (with large (5/3) . This would explain why these parameter values lead 
to significantly increased dissipation. Further studies of this problem are clearly 
needed if we are to understand the r-modes in a superfluid star fully. Investiga- 
tions of this problem should be strongly encouraged, both in Newtonian theory 
and general relativityllj. 

4.5. Ekman layer in neutron stars with a crust 

In addition to the core superfluid, we need to consider effects due to the 
presence of a solid crust in an old neutron star. The melting temperature of 
the crust is usually estimated to be of the order of lO^'' K (for a non-accmting 
neutron star), so the crust may form shortly after the neutron stars birthEj. 

That the presence of a solid crust will have a crucial effect on the r-mode mo- 
tion can be understood as follows: Based on the perfect fluid mode-calculations 
we anticipate the transverse motion associated with the mode at the crust-core 
boundary to be large. However, if the crust is assumed to be rigid the fluid mo- 
tion must essentially fall off to zero at the base of the crust in order to satisfy 
a no-slip condition (in the rotating frame of reference). We can estimate the 
relevance of the crust using viscous boundary layer theory0. The region im- 
mediately beneath the crust then corresponds to a so-called Ekman layer. The 
thickness of the boundary layer {5) can be deduced by balancing the Coriolis 
force and shear viscosityc^l: 



where rj is the shear viscosity coefficient (jT^)- As is easy to show, this would 
correspond to a few centimetres for typical parameters of a rapidly rotating 
neutron star. In other words, 5 << R and the boundary layer approximation is 
warranted. 

The dissipation timescale due to the presence of the Ekman layer roughly 
follows from 

tEk - (94) 

where Re — pbR^Vl/rj is the Reynolds number (the ratio between the Coriolis 
force and viscosity), and i?b and are the location of, and density in, the Ekman 
layer, respectively. For typical neutron star parameters one would expect the 
base of the crust to lie at pt, « 1.5 x 10^"* g/cm . We would then have Re ~ 10^^, 
and we see that the presence of a solid crust leads to a dissipation channel that 
is many orders of magnitude more effective than the standard shear viscosity. 

A more detailed estimate of this effect has been made by Bildsten and 
UshomirskyEJ. The required dissipation rate follows from an integral over the 
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surface area at the crust-core boundary 

|(5u|2 



dE 
'dt 



Ek 



(95) 



For canonical neutron star parameters and an n = 1 polytrope we can estimate 
that the base of the crust corresponds to Rb ~ 9.4 km. When combined with 
the estimated mode-energy, this leads to an estimated damping timescale (note 
that this estimate in valid only for a i? = 10 km and M = IAAIq star) 



tEk » 1.4 X WTqP_ 



1/2 



(96) 



1.0-- 
0.8- 

0.6- 

Oh 

0.4- 
0.2- 



0.0-1 ■ , ■ , ■ , ■ , ■ , ■ 1 

5 6 7 8 9 10 11 

l0g,„ (T/IK) 

Fig. 6. The r-mode instability window for a canonical fluid neutron star (lower solid curve) 
is compared to that for a star with an Ekman layer operating at the base of the solid crust 
(upper solid curve). Also shown (as dashed curves) are the corresponding instability windows for 
a constant density star. It should be noted that the strength of the r-mode instability is only 
weakly dependent on the stiffness of the equation of state. 

This estimate was arrived at by appealing to an analogy between the problem 
at hand and the standard one of an oscillating plate in a viscous fluidEil. One 
would not expect the change in geometry and inclusion of the Coriolis force in 
the r-mode problem to change the order of magnitude of the effect. That this 
assumption is correct has been shown by RieutordS, who used an exact solution 
to the incompressible problem to show that when the angular dependence of the 
velocity in the Ekman layer is accounted for the dissipation timescale for the 
I — m = 2 r-mode is a factor of 1.74 shorter than (|9^). The resultant estimate 
of the dissipation due to the presence of the Ekman layer is 

<Ek ~ 830TgP%^ s . (97) 

We compare the r-mode instability window obtained from this estimate to the 
normal fluid result in Figure 0. Clearly, a solid crust has an influence that is 
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far greater than any previously considered dissipation mechanism. For example, 
our estimate suggests that all neutron stars with a crust are stable at rotation 
periods longer than roughly 5 ms. 

The crust-core interface has been the focus of several recent studies. These 
studies add further dimensions to the problem. The interplay between the r- 
modes in the fluid core and modes in the solid crust is particularly interesting. 
It has long been known that the crust supports toroidal modes of oscillation 
(whose frequencies depend directiy on the crust shear modulus )EZll3. Detailed 
calculations by Yoshida and Leel£3 (see also Levin and UshomirskyE3) show that 
as the spin of the star increases the unstable r-modes will undergo a series of 
so-called avoided crossings with the crust modes. In a nonrotating star, the 
fundamental I = 2 crust mode has a frequency 



u;,{n^Q)^m-^^^ . (98) 

In the rotating case the mode-frequency changes in such a way that 

777 

co^in) « co^in = 0) + (99) 

(in the rotating frame). Now it is easy to see that the toroidaLpode becomes 
comparable to the r-mode frequency at a rotation frequencyll3H 

f^cross » 5 X 10-2^!^ (100) 

Hence, avoided crossings between the modes occur at very low rates of rotation. 
This is conceptually interesting and it may have repercussions on many of our 
estimates regarding the r-mode instability. Most importantly, it casts doubt on 
extrapolations of the slow-rotation results into the regime of fast rotation. 

In addition to studying this new feature of the problem. Levin and UshomirskyE3 
show that the assumption of a rigid crust, which was made in the above esti- 
mates of the Ekman layer dissipation rate, is likely not warranted. They show 
that the r-mode typically extends into the crust. We can quantify the extent 
to which this affects the estimated dissipation timescale ( p7| ) in terms of the 
ratio A = |At;|/|u|, where Av is the difference between the core fluid velocity 
immediately beneath the crust and the velocity induced in the crust and v is the 
fluid velocity in absence of the crust. The rigid crust assumption corresponds 
to A = 1. In contrast, the toy model calculations of Levin and UshomirskyEj 
suggest that a typical results should be A « 0.1. This then affects the Ekman 
layer timescale by a factor of 1/A^, so the true dissipation may well be a factor 
of at least 100 weaker than (p7|). 

Despite some recent advances in our understanding of the effects that a solid 
crust may have on the r-modes it is clear that several crucial issues remain to 
be investigated in detail. For example, the inner crust of a neutron star (out to 
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the neutron drip density) will likely be permeated by superfluid neutrons. It is 
not at all clear at the moment whether one should expect these neutrons to be 
strongly pinned to the crust nuclei or not. But if the superfluid is at all free to 
move relative to the crust it will likely lead to a weaker Ekman layer dissipation 
on the r-modes. 

4.6. Strange/quark stars 

Ever since Wittenll suggested that strange matter (a conglomerate of strange, 
up and down quarks) would be the most stable form of matter, it has been specu- 
lated that strange stars might exist in the universe0il'H It has been suggested 
that some (perhaps all) pulsars are actually strange stars. The observational 
evidence for this is, however, rather tenuous. The main reason for this is that 
strange stars are held together by both the strong interaction and gravity, and 
at a mass of one and a half solar masses, gravity dominates over the strong 
interaction. Thus it is not easy to distinguish a strange star from a neutron 
star, particularly iiot since the strange star is expected to be "covered" by a 
thin nuclear crustEH. Still, the existence of strange stars remains a conjecture 
which is consistent with, for example, millisecond pulsar observations. 

As was pointed out by MadsenE2l, the r-mode instability may provide the 
means for distinguishing between strange stars and neutron stars. The main 
reason for this is that the viscosity coefficients are rather different in the two 
cases. While the shear viscosity in a strange star would be comparable to that of 
a neutron star, the bulk viscosity would be many orders of magnitude stronger 
than its neutron star counterpart. This has interesting effects on the r-mode 
instability. 

For the shear viscosity the relevant coefficient would be 

77 « 1.7 X 10i« (^^) plt^X^'/'g/cms (101) 
where is the fine-structure constant for the strong interaction. This leads to 
7.4 X 10^ (^) M^j/'Rll^'T^/' s (102) 

Meanwhile, the situation is more complicated for the bulk viscosity (which 
now is a result of the change in concentration of down and strange quarks in 
responaato the mode oscillation). The relevant viscosity coefficient takes the 
fornffl 

where the coefficients a and /3 are given by Madsen. From this we can imme- 
diately deduce that the bulk viscosity becomes less important at both very low 
and very high temperatures. For low temperatures we find that 

47 « 7.QMl,R^,^Pl,T,''m^,% s (104) 
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where rriioo represents the mass of the strange quark in units of 100 MeV. Mean- 
while, at high temperatures we cannot readily write down an expression with 
the appropriate dependence on M, R, P etcetera. Instead, we have to perform 
numerical calculations for each stellar model. Typical results (for canonical pa- 
rameters) are shown in Figure ^. From this figure we immediately see that the 
r-mode instability would not - contrary to the case for neutron stars - be active 
in strange stars with a core temperature of 10^ K. In a strange star the r-modes 
are unstable at lower temperatures (between 10^ — 5 x 10^ K) and also at tem- 
peratures above a few times 10^ K. This is interesting since it means that the 
instability would be active for a brief period after a strange star is born. Then 
the mode would become stable until thousands of years later when the star has 
cooled sufficiently to enter the low-temperature instability window. 
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Fig. 7. The r-mode instability window for a strange star (solid line). As comparison we show the 
corresponding instability results for normal fluid and crusted neutron stars (dashed curves) from 
Figure ^ The decreased importance of bulk viscosity in hot strange stars is notable. 

Similar estimates_would be valkLfor a neutron star with a quark core, a 
so-called hybrid staii23 and MadsenEll has commented on the viscosity of some 
rather exotic states of matter that may be relevant for hybrid star cores. 

5. Is the r-mode instability astrophysically relevant? 

In the previous sections we have described the nature of the r-modes of 
a rotating neutron star. We have discussed the associated gravitational-wave 
driven instability and shown that it may operate in a neutron star with core 
temperature in a window between roughly lO'^ — 10^° K. This means that the 
instability should be relevant for hot young neutron stars that are born spinning 
at a rate above perhaps 5% of the dynamical Kepler limit. 

We now want to discuss the possible astrophysical consequences of the r- 
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mode instability. By necessity this means that we will extrapolate the available 
results considerably. Astrophysical neutron stars are much more complex than 
the simple models for which the properties of the r-modes have so far been 
studied. Nevertheless, the speculations that we describe in this section are 
interesting and potentially of great importance. After all, they may eventually 
turn out to be largely correct! We must, of course, keep in mind that our 
understanding of much of the involved physics is (at best) rudimentary. 

5.1. A phenomenological spin-evolution model 

If we want to make the discussion quantitative rather than qualitative we 
must model the actual evolution of an r-mode as it grows and spins down the 
star. Given the current models, this is a challenging task and any model we 
devise is likely to be questionable. Still, it is important that we try to gain 
insight into this issue. After all, once we have an initial prediction we can try 
to identify the crucial elements of the model and work to improve on them. It 
is also relevant to address many important astrophysical questions at an "order 
of magnitude" level. 

While the early growth phase of an unstable mode can be described by linear 
theory, an understanding of many effects (such as coupling between different 
modes) that will become important and may eventually dominate the dynamics 
require a nonlinear calculation. While such studies are outstanding we can only 
try to capture the essential features of the behaviour. Intuitively one would 
expect a growth of an unstable mode to be halted at some amplitude. As the 
mode saturates it seems plausible that the excess angular momentum will be 
radiated away and the star will spin down. This general pictur e is supported 
by detailed studies of instabilities in rotating ellipsoids0'@O@. 

We model the evolution of a spinning star governed by the r-mode instability 
using the phenomenological two-parameter model devised by Owen et al.lll That 
is, we consider the spin rate Q and the mode-amplitude a as our key quantities. 
Recall that the latter is defined by expressing the velocity field (measured by 
an observer that is co-rotating with the star) as (for I = m modes) 



Assume that the total angular momentum of the system can be decomposed 

as 



where the first term represents the bulk rotation of the star and the second 
is the canonical angular momentum of the r-mode. This representation makes 
sense since Jc corresponds to the second order change in angular momentum ^ 

ITo be more precise: There is an ambiguity in deciding whether Jc is the second-order change 
in angular momentum, because that second-order change also includes an arbitrary second 




(105) 



J = m + Jc 



(106) 
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due to the presence of the modalj'ESl (m absence of viscosity or radiation) . In 
the case of r-modes we have 



Jc 



Sfla^JMR'^ 



where we have defined 



.1 f'^ 6,j 3/287r n^Q 
" An¥ Jo - \ 1.635 X 10-2 n = 1 



(107) 



(108) 



We now assume that angular momentum is only dissipated via gravitational 
waves. The gravitational-wave luminosity follows from the I — m = 2 current 
multipole formula, and we get 



JMR-^ 

tavj 



Taking a time-derivative of (106) and defining 

2/5 : 

0.261 n=l 



we have 



1 a J — iilaJ— = ^ 

2 dt dt MR^ 



(109) 



(110) 



(111) 



Here we have neglected the dependence of the moment of inertia on the spin 
rate (terms of order fi^). It is easy to include such terms and show that they 
have a minor effect. 

A second equation follows from the expression for the energy of the mode 



(as measured in the rotating system), cf. (58) 

1 



E 



a^n^MR^J 



(112) 



and the fact that the mode grows according to 



dE 

— = -2E , 
dt \ t 



1 



where 



1 11 , . . 

— \ h other dissipation terms 

^diss ^sv ^bv 



(113) 



(114) 



order addition of an solution to the time-independent linearized equations that simply adds 
differential rotation to the equilibrium star. This ambiguity can be resolved by requiring that 
the »ie]|[,U|L bation be canonical and that the Lagrangian displacement have no second-order 
part eStEj. Under these conditions Jc is the second order change in the angular momentum. 
For a growing mode, requiring that the second-order part grows exponentially in the manner 
e't/r), (with t the growth time of the mode energy) eliminates the ambiguity in a physically 
more appropriate way, but the consistent solution requires one to include the also the second 
order radiation-reaction term. 
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These two equations combine to give 



dt dt \ tgw ^diss , 



We can now combine (111) and (115) to get "evolution equations" for the 



mode amphtude and the bulk rotation rate: 

da a A_ 1 

— = a- (116) 

ISS 

and 

§ = -^3a'nj-l- (117) 

at A+ tdiss 

where 

A±=i±^a^J (118) 

These equations govern the mode-evolution in the phase when the amplitude 
of the r-mode grows, as well as in the late phase when the temperature has 
decreased sufficiently to again make the mode stable. 

Since they should be adequately described by linear perturbation theory, the 
very early and late parts of the evolution of an unstable mode are comparatively 
well understood. Unfortunately, it is during the nonlinear phase that the main 
spin-down is anticipated to occur. Hence, we need to model the nonlinear regime 
in a "meaningful" way. As our starting point we adopt the simple view that the 
mode- amplitude will saturate at a critical value as < 1. The assumed upper 
limit of the saturation amplitude, as — 1 corresponds to the r-mode carrying 
a considerable part of the angular momentum of the system. We can readily 
estimate that 



Jc 



O.la^ (119) 



21 

Hence, a value of as « 1 should be considered as very large (cf. also (^)). 
We assume that a stays constant throughout the saturated phase and that fl 



evolves according to ( |lll|) with a = 0. This means that we have 

dfl 1 arja^j 



dt A— ^gw 



(120) 



Since the viscosities are sensitive to temperature changes we need to model 
also the thermal evolution of the star in some suitably simple way. To do 
this we account for two main processes: i) cooling due to (say) the modified 
URCA process, and ii) reheating due to energy deposited from the mode to 
the heatbath via the shear viscosity. In stars with a solid crust we also include 
energy dissipated into heat in the Ekman layer. The bulk viscosity leads to the 
creation of neutrinos that can safely be assumed to escape from the star on a 
very short time-scale at the temperatures we consider. 
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We will typically initiate the evolution of the star at some temperature (recall 
that neutron stars are born above 10^^ K and cool to 10^ K in a few months). 
To relate the initial temperature to the thermal energy available in the heat 
bath we can useEB 

Et = 3.5 X 10^'^ MiAPi5^^T^ erg (121) 

Then we model the subsequent thermal evolution of the star by assuming that 
Et evolves according to 

— TT" = '^^ly + -E-visc (122) 

at 

where the luminosity due to the modified URCA process ifl 

^VRCA ^ ^ ^ ^ 104OMi.4p}f Tg" erg/s (123) 

This mechanism is anticipated to dominate the cooling in the relevant range of 
temperatures, but one should not forget the jipssibility of more rapid cooling 
due to for example the direct URCA reactions^. The possible outcomes should 
a faster cooling mechanism be operating_have been discussed by Andersson, 
Kokkotas and SchutzcJ and Yoshida et al.E!l 

The reheating due to shear viscosity is readily determined from 

2an^MR^J , , 

Es. = 7 (124) 

^sv 

The heat generated due to the presence of an Ekman layer in an older neutron 
star with a crust follows from a similar expression. 

5.2. The spin-evolution of young neutron stars 

Given the general features of the r-niode instability we expect that it could 
be relevant in newly born neutron stars. Hence it is interesting to consider the 
observational data regarding the initial spin rate of young pulsars. The best 
studied young pulsar is the Crab (PSR053H-21), whose initial period is esti- 
mated_( assuming the standard magnetic braking model) to have been about 
19 msLj. Even the recently discovered young 16 ms X-ray pulsar in the sup er- 
nova remnant N157B probably had an initial period no shorter than a few m^ 
(assuming a braking index typical of young pulsars). These estimates are in 
clear contrast with the shortest known period of a recycled pulsar of 1.5fi ms, 
and with the theoretical lower limit on the period of about 0.5 to 2 mai3, de- 
pending on the equation of state. The available data essentially suggests that 
neutron stars are formed spinning slowly, at perhaps Q/^Ik < 0.1. From a 
theoretical point of view, this is rather surprising since, assuming that angular 
momentum is conserved in the collapse event that forms the neutron star, young 
pulsars ought to be spinning close to the Kepler limit. 

One possible explanation for the slow i niti al spin rate of newly born neutron 
stars was proposed by Spruit and Phinney0, who argue that magnetic locking 
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between the core and the envelope of the progenitor star may prevent the col- 
lapsing core from spinning rapidly. Their estimates suggest that the core would 
actually rotate far too slowly to lead to (by conservation of angular momen- 
tum) young neutron stars spinning as fast as they do. To explain the observed 
rotation rates, Spruit and Phinney propose that the neutron star spin is due 
to the birth kicks that also (or alternatively) produce the large linear velocities 
observed for pulsars. 

Even if correct, Spruit and Phinney 's model still suggests that some pulsars 
are born spinning rapidly. Theoretically, it is easy to estimate that birth kicks 
that lead to velocities larger than 500 km/s can also produce rotation at (or 
above) the Kepler limit. Furthermore, the 16 ms pulsar in N157B is clear 
evidence that some pulsars form with periods shorter than say 10 ms. For 
such young neutron stars, the secular r-mode instability may play a role in 
determining the rotation rate. It is this possibility that we focus on here. 

As soon as we bring the spin-evolution model of the previous section to bear 
on the problem we see that the key parameter is the saturation amplitude a^. 
Provided that as is sufficiently large, the r-modes will spin a young neutron 
star down appreciably. If we assume that the neutron star is born with core 
temperature well above 10^" K and that it initially spins at the Kepler limit, 
the r-mode instability comes into play within a few seconds as the star cools 
and enters the instability window. The mode then grows from some small 
initial amplitude to the saturation level in a few minutes. Once the mode has 
saturated, the star spins down. At some point the star has cooled (or spun 
down) sufficiently that the r-mode is again stable. Then the mode amplitude 
decays and the star presumably enters a phase where magnetic braking takes 
over and dominates the spin-evolution. Examples of this scenario are shown in 
Figure ||. 

The final spin period depends on several factors. The most important of 
these are the saturation amplitude, the cooling rate and whether a crust forms 
during the evolution. If we consider a simple perfect fluid model the r-mode 
spin-down leads to P ~ 12 — 22 ms for as in the range 0.01 — 1 (and canonical 
neutron star parameters) . This result is clearly consistent with observations for 
many young pulsars (in particular the Crab). 

As we discussed in Section 4.5, the r-mode instability may be strongly sup- 
pressed in neutron stars with a solid crust. Given that the melting temperature 
of the crust could be as high as 10^° K, the crust may form shortly after the 
neutron star is born and we need to discuss the effect that this may have on the 
r-mode scenario. The interplay between a growing mode and the formation of 
the crust leads to complicated questions that require much further study. For 
example, it is not at all clear to what extent the melting temperature of the 
crust is affected by large scale surface waves in the star. Here we will not address 
such questions in detail (but see comments in Section 6). Instead we consider 
what may be a "worst case scenario" as far as the instability is concerned. We 
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Fig. 8. Period evolutions resulting from the r-mode spin-down scenario. We assume that the 
star is isothermal and that it cools entirely due to the modified URCA process. Two dilTerent 
evolutions corresponding to an r-mode saturated at = 0.01 and 1, respectively, are shown as 
solid lines. The dashed curve indicates the evolution that would follow if the r-mode was initially 
excited to the saturation level rather than given a very small initial amplitude (a = 10~^) at the 
onset of the instability. It is notable the initial amplitude has little effect on the period at the end 
of the spin-down phase. 

assume that the r-mode will not be able to prevent the crust from forming (and 
that the melting temperature is at the high end of the anticipated range) . The 
main dissipation mechanism is then due to the presence of the Ekman layer (at 
least above the superfluid transition temperature). Resulting spin evolutions 
are shown in Figure |[ In this scenario, the period reached after the spin-down 
phase lies in the range P « 2.5 — 4.5 ms. Just as in the perfect fluid case, this is 
an interesting prediction since it agrees quite well with the data for the recently 
discovered 16 ms pulsar PRS J0537-6910 in N157B. 

We have thus seen that the r-mode instability can lead to young neutron 
stars being spun down to rotation rates that would agree quite well with ex- 
trapolations from current observations. Furthermore, we have predicted that 



scenario applies depends sensitively on the early cooling of the star and the 
crustal formation temperature. To illustrate this, we consider the evolution of 
a neutron star just after its birth in a supernova explosion. We might expect 
to model its r-mode evolution using the normal (crust-free) fluid viscous damp- 
ing times for stellar temperatures above the melting temperature of the crust 
(Tm), and the viscous boundary layer damping time for temperatures below Tm- 
However, the situation is probably a little more complicated than this. Recall 
that the latent heat (i.e. the Coulomb binding energy) of a typical crust is 



hot young neutron stars may follow two different 
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£;iat»10*« erg 



(125) 
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Fig. 9. Spin-down evolutions for a neutron star in which a solid crust forms before the r-mode 
has grown appreciably. The presence of an Ekman layer at the base of the crust leads to a strong 
dissipation of the mode energy and the final spin-period after the instability phase is much shorter 
than in the crustless case. 



As is clear from (58) the energy in the r-mode (which grows exponentially on 
a timescale ig^ ~ 20 — 50 s) will easily exceed i?iat, provided that the time 
taken for the star to cool to is sufficently long. This would then prevent the 
formation of the crust, even when T < Tm- Then the star will spin down in the 
manner shown in Figure ||. This phase will end either because the star leaves 
the instability window or because the mode energy has fallen below the crustal 
binding energy. Supposing that the r-mode energy is distributed in such a way 
that (say) 10% is located in the region where the crust will form (the outer km 
or so of the star), we reach balance, i.e. i^iat ~ O.liJmodo, when 



10a, 



(126) 



In other words, for a saturation amplitude of order unity the crust would be 
unable to form until the that has spun down to a rotation rate of roughly 100 Hz. 
Again, this would be consistent with the extrapolated initial spin rates of many 
young pulsars. On the other hand, if the mode is not given time to grow very 
large it will not prevent crust formation at Tm- The the r-mode instability 
would spin the star down as in Figure ^. A likely key parameter (in addition 
to the cooling rate, the melting temperature and the extent to which a large 
amplitude mode can melt the crust, see section 6) is whether the supernova 
collapse leads to a large initial r-mode amplitude or not. 

Besides explaining the slow rotation of young pulsars, the r-mode instability 
has the consequence that young neutron stars can only reach rotation periods 
shorter than (say) 3-5 ms if they are recycled by accretion in a binary system. 
The alternative model, that these stars are formed by accretion-induced collapse 
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of a white dwarf, is inconsistent with the r-mode scenarioEa. The cohapse would 
form a star hot enough to spin down because of the instabihtjf However, the 
situation would be quite different in the case of strange starsL^. As shown in 
Figure ^, the instability window for strange matter differs significantly from the 
neutron star one. This means that the r-mode instability would be less effective 
(apart from at very high temperatures) for strange stars. Furthermore, a strange 
star cools much slower than a neutron star. Assuming that the strange star cools 
on a timescale t ~ lO^^Tg"* yrs, we find that the strange star reaches P « 3 ms 
after 4 x 10^ yrs. This means that an observation of a young millisecond pulsar 
spinning faster than (say) 3 ms may indicate the presence of a strange staiO. 

5.3. Are the gravitational waves observable? 

Having suggested that the r-mode instability may spin a newly born neutron 
star down to a fraction of its initial spin rate in a few months we want to know 
whether the gravitational waves that carry away most of the stars initial angular 
momentum are detectable. This is a particularly relevant question given the 
generation of large-scale interferometers that is about to come into operatiorfl. 

We assess the detectability of the emerging gravitational waves in the stan- 
dard way. First of all, we note that the frequency of the emerging gravitational 
waves is (for the main I = ni = 2 r-mode) 

/gw = ^ . (127) 

From the gravitational- wave luminosity we readily deduce that the dimensionless 
strain amplitude follows from 

h{t) = 7.54 X 10-23aJAfi.4i??o (^^^^^ (128) 

where D is the distance to the source, here assumed to be in the Virgo cluster. 
At this distance one would expect to see several neutron stars being born per 
year. The typical r-mode scenarios shown in Figures || and ^ then lead to the 
gravitational-wave strains illustrated in Figure ^ Clearly, this amplitude is 
not sufhciently strong that it can be observed without a detailed data analysis 
strategy. To assess the possible improvement that would follow if such a strategy 
could be developed for the r-mode signal we will use the standard matched- 
filtering approach. This gives us an idea of the improvement that a specially 
tailored data analysis approach may bring to the r-mode detection problem, even 
though it must be rec ogn ized that matched filtering is unlikely to be possible 
for this kind of signals I^IhI. 



With the Fourier transform of the gravitational-wave signal defined by 

/C30 
h{t)e^''f'dt (129) 
-OO 
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Fig. 10. The dimenslonless gravitational-wave amplitude following from a typical r-mode sce- 
nario. The two phases of evolution (the early growth phase of the mode and the spin-down after 
saturation) are easily distinguished. The second phase is the most likely to be observed with future 
detectors. We show signals corresponding to both fluid neutron stars (thick solid line) and ones 
with a crust (dashed). In both cases we assume that the saturation amplitude is as = 1. It should 
be noted that the signal lasts 10^ s in the fluid case but only 10^ s for the crusted star. 



the signal-to-noise ratio attainable by optimal filtering can be determined from 



rrns 



(130) 



where hrms = \/ fShif), Sh being the spectral density of the strain noise in the 
detector. Now introducing the characteristic amplitude he we can write 



^1 ^2 

N 



f 



where 



(131) 



(132) 



The last relation follows via the stationary phase approximation. It can be inter- 
preted as meaning that the detectability of an almost periodic signal improves 
as the square root of the number of cycles radiated in the time it takes the fre- 
quency to change by /. In Figure |l^ we compare this characteristic amplitude to 
estimated sensitivity curves for the new generation of interferometric detectors. 
We immediately see that the r-mode instability would be a promising source 
for an advanced generation of detectors such as LIGO II. One can estimate that 
the gravitational waves fro m a hot young neutron star will be detectable with a 



good signal-to-noise rati' 



000 



lo; 



This makes the r-mode instability one of the 
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most promising sources of detectable gravitational waves. It is unlikely that the 
r-modes will be observed by the first generation detectors, however. For this to 
happen requires a unique event in our Galaxy or the local group. One would 
typically expect a supernova in the galaxy every 30 years or so, which means 
that we would be extremely lucky to see such an event. Anyway, the chances of 
future detection are promising and they may improve as our understanding of 
the instability becomes more detailed and better theoretical templates for the 
gravitational-wave signal can be constructed. 




f(Hz) 



Fig. 11. The effective amplitude achievable after matched filtering (which is somewhat unre- 
alistic for the r-mode instability) compared to the expected sensitivity of the new generation of 
interferometric detectors. The two gravitational-wave signals are the same as in Figure id. 



In addition to signals from individual neutron stars, there will be a stochas- 
tic background of gravitational waves from cosmological neutron stars spinning 
down due to the instability. Estimates show that this background will be diffi- 
cult to detect, but it could perhaps be relevant for an advanced generation of 
detectorSH. 



6. Accreting neutron stars 

In the last few years the Rossi X-ray Timing Explorer has provided a wealth 
of observational data regarding accreting neutron stars. These observations 
present a challenge for theorists in that neutron stars in low-mass X-ray bina- 
ries (LMXB) seem to be confined to a rather narrow range of spin frequencies, 
perhaps 260-590 HzEHj. Different models have been proposed to explain this 
surprising result. The first model (due to White and ZhangE£3) is based on the 
standard magnetosphere model for accretion induced spin-up, while two other 
models are inspired by the idea that gravitational radiation balances th e ac cre- 



tion torque. In the first such model for the LMXB (proposed by BildstenEill and 
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recently refined by Ushomirsky, Cutler and Bildstenll23), the gravitational waves 
are due to an accretion induced quadrupole deformation in the deep neutron 
star crust. The second gravitational-wave model relies on the r-modejnstability 
to dissipate the accreted angular momentum from the neutron start2j. 

The idea that gravitational waves from unstable r-modes provide the agent 

torque was first analyspd in detail by Andersson, 



that balances the accret 



on 



Kokkotas and Stergioulasll23 (but see also Bildstenll2Zl) . Originally, it was thought 



that an accreting star in which the r-modes were excited to a significant level 
would reach a spin-equilihrium, very mucL in the vein of early suggestions by 



Papaloizou and Pringleti^l and Wagoner — . Should this happen, the neutron 
stars in LMXB would be prime sources for detectable gravitational waves. That 
the associated (essentially periodic) gravitational waves would be detectable can 
be seen in the following way: Assume that the r-mode instability is active and 
provides a limit on the spin of Sco XI, which is the strongest X-ray source in the 
sky and therefore a prime candidate for this kind of speculation. If we assume 
that the average accretion rate onto the neutron star is M « 3 x lO~^M0/yr 
and that the accretion torque is balanced by gravitational radiation we deduce 
that h « ,3. 5 X 1fl~^^ he ^ 10^^^ after two weeks worth of signal has been 
accumulatedll23. 

However, this idea is probably not viablel* In addition to generating gravi- 
tational waves that dissipate angular momentum from the system, the r-modes 
will heat the star up (via the shear viscosity and the Ekman layer). Since the 
viscosity gets weaker as the temperature increases, the mode-heating triggers a 



thermal runawaylliSlli^ and in a few months the r-mode would spin an accreting 
neutron star down to a rather low rotation rate. , This essentially rules out the 
r-modes in galactic LMXB as a source of detectable gravitational waves, since 
they will only radiate for a tiny fraction of the systems lifetime. 



This mechanism could still be of great astrophysical significanceE£ll. Hence, 
we want to model how the potential presence of an unstable r-mode affects the 
spin-evolution of rapidly spinning, accreting neutron stars. To do this we use 
the phenomenological two-parameter model described in section 5.1. Then the 
following picture emerges: After accreting and spinning up for something like 
10^ years, the star reaches the period at which the r-mode instability sets in. 
For a canonical neutron star this corresponds to a period of 1.5 ms (at a core 
temperature of 10^ K). It is notable (albeit likely coincidental) that this value 
is close to the 1.56 ms period of the fastest known pulsar PSR1937-f 21. Once 
the r-mode becomes unstable (point A in Figure [T^ ), viscous heating (mainly 
due to the energy released in the Ekman layer) rapidly heats the star up to a 
few times 10^ K. The r-mode amplitude increases until it reaches the saturation 
level (amplitude as), at which unspecified nonlinear effects halt further growth 



§ Although it should be noted that the Wagoner scenario might be relevant if accretion at 
extreme rates is possible. The role of the unatetHle r-modes during such hypercritical accretion 
has been discussed by Yoshida and Eriguchiti3. 
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Fig. 12. The r-mode instability window relevant for old neutron stars. We show results for the 
simplest (crust-free) model (thin solid line), as well as for a star with a crust (thick solid line). We 
illustrate two typical r-mode cycles (for mode saturation amplitudes ag = 0.1 and 1), resulting 
from thermo-gravitational runaway after the onset of instability. Once accretion has spun the star 
up to the critical period (along the indicated spin-up line (thick vertical line)) the r-mode becomes 
unstable and the star evolves along the path A-B-C-D. After a month or so, the mode is stable and 
the star will cool down until it again reaches the spin-up line. For comparison with observational 
data, we indicate the possible range of spin-periods inferred from current LMXB data (shaded 
box) as well as the observed periods and estimated upper limits of the temperature of some of the 
most rapidly spinning millisecond pulsars (short arrows). 



(point B in Figure |T^). Once the mode has saturated, the neutron star rapidly 
spins down. When the star has spun down to the point where the mode again 
becomes stable (point C in Figure ^2|), the amplitude starts to decay and the 
mode eventually plays no further role in the spin evolution of the star (point D 
in Figure Two examples of such r-mode cycles (corresponding to as = 0.1 
and 1, respectively) are shown in Figure [l^. 

This simple model leads to some interesting quantitative predictions. First 
of all, the model suggests that an accreting star will not spin up beyond 1.5 ms. 
This value obviously depends on the chosen stellar model, but it is independent 
of the r-mode saturation amplitude and only weakly dependent on the accretion 
rate (through a slight change in core temperature). In fact, the accretion rate 
only affects the time it takes the star to complete one full r-mode cycle. As soon 
as the mode becomes unstable the spin-evolution is dominated by gravitational 
radiation and viscous heating. Once the star has gone through the brief phase 
when the r-mode is active it has spun down to a period in the range 3-5 ms 
(corresponding to 0.01 < < 1). At this point the mode is again stable 
and continued accretion may resume to spin the star up. Since the star must 
accrete roughly O.IMq to reach the instability point, and the LMXB companions 
typically have masses in the range 0.1 — OAMq, it can pass through several "r- 
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mode cycles" during its lifetime. 

We can readily confront these results with current observations^. We note 
that the model predicts that an accreting neutron star must, once it has been 
spun up beyond (say) 5 ms, remain in the rather narrow range of periods 1.5 — 
5 ms until it has stopped accreting and magnetic braking slows it down. Since 
a given star can go through several r-mode cycles before accretion is halted 
one would expect most neutron stars in LMXB and the millisecond pulsars to 
be found within this range. As is clear from Figure this prediction agrees 
well with the range of rotation periods inferred from observed kHz quasiperiodic 
oscillations in LMXB. The observed range shown in Figure ^2] corresponds to 
rotation frequencies in the range 260- 590 Hzlii!3. The model also agrees with the 



observed data for millisecond pulsarsliiliJ, which are mainly found in the range 



1.56 — 6 ms. In other words, the r-mode runaway model is in agreement with 
current observed data for rapidly rotating neutron stars. 

However, many crucial pieces of physics have not yet been included in this 
model. For example, we based the discussion on the simplest estimate (|97| ) of the 
Ekman layer dissipation, rather than morerefined ones attempting to include 
the crust-core interaction in a detailed wayO'EJ. Furthermore, we assumed that 
the heat released in the Ekman layer will not be able to melt the crust. This 
latter issue ha. s recently been studied in some detai l by Lindblom, Owen and 
Ushomirskyllia (see also Wu, Matzner and ArraaliJ) , who solve the appropriate 
heat equation incorporating thermal conductivity for the core fluid and the crust 
as well as neutrino cooling due to the modified URCA process. They find that, 
for a star spinning at the Kepler limit, an r-mode with an amplitude exceeding 
a ~ 5 X 10^'^ will be able to melt the crust. This means that if the r-mode 
instability is active in a young neutron star, the formation of the crust may 
be significantly delayed (as we argued in section 5.2). But it is not yet clear 
what will happen if the r- modes do indeed manage to melt an existing crust. 
As soon as the crust melts the Ekman layer (that led to the excessive heating) 
disappears and the material will rapidly cool down to ajevel where the crust 



would begin to form again. As argued by Lindblom et alt-Lj, the likely outcome 



is a solid-fluid mixed state which will be very difficult to model in detail. 

We conclude this section by discussing briefiy the detectability of the gravita- 
tional waves that are radiated during the relatively short time when the r-mode 
is saturated and the star spins down. Since the r-mode is active only for a 
small fraction of the lifetime of the system (something like 1 month out of the 
10^ years it takes to complete one full cycle) the event rate is far too low to make 
galactic sources relevant. However, it interesting to note that the spin-evolution 
is rather similar to that of a hot young neutron star once the r-mode has reached 
its saturation amplitude. This means that we can analyse the detectability of 
the emerging gravitational waves using the framework of section 5.3. We then 



'Here it is worth mentioning that Brown and Ushomirsky|li£ 



have discussed how observations 



of X-ray transients may provide a constraint on the r-mode amphtude in such systems. 
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find that these events can be observed from rather distant galaxies. For a source 
in the Virgo cluster these gravitational waves could be detected with a signal to 
noise ratio of a few using LIGO II. However, even at the distance of the Virgo 
cluster these events would be quite rare. By combining a birth rate for LMXB 
of 7 X 10~^ per year per galaxy with the fact that the the volume of space 
out to the Virgo cluster contains ~ 10'^ galaxies, and the possibility that each 
LMXB passes through (say) four r-mode cycles during its lifetime we deduce 
that one can only hope to see a few events per century in Virgo. In order to see 
several events per year the detector must be sensitive enough to detect these 
gravitational waves from (say) 150 Mpc. This would require a more advanced 
detector configuration. 

7. Steps towards more detailed modelling 

From the discussion in the previous sections it should be clear that the 
present results indicate that the unstable r-modes may have great astrophysi- 
cal significance, and that the associated gravitational-wave signals may well be 
detected within a few years. These ideas are obviously very exciting, but it is 
important to realize that our modelling must be improved in many directions 
before any reliable conclusions can be drawn. In this section we discuss some 
recent work that provide extensions (or alternatives) to the ideas that we have 
already described. 



7.1. Inertial modes and rapid rotation 

As was mentioned in the introduction, the r-niodes belong to a larger class 
of pulsation modes which in fluid mechanics are commonly described as "inertial 
modes"E3. These modes are all dominated by the Coriolis force and for slowly 
rotating stars their frequencies scale linearly with fi. In non-isentropic stars the 
rotation dominated modes are all toroidal in nature (the r-modes). In contrast, 
the space of zero frequency modes of isentropic spherical stars includes the 
spheroidal y-modes as well as the r-modes. This large degenerate subspace 
of zero- frequency modes is then split by rotation, and the inertial modes of 
isentropic stars generally have a "hybrid" natureEj (their spherical limit is a 
mixture of toroidal and spheroidal perturbations). 

Detailed calculations of such hybrid modes in Newtonian stars h ave been 



made by Lockitch and FriedmanE^I as well as Yoshida and LeeE£llii3 (see also 
Lee et alE£2l and the studies of Maclaurin spheroids discussed below). These 
studies show that the low frequency mode spectrum of a rotating star is tremen- 
dously rich. Furthermore, many of the hybrid modes are affected by the CFS 
instability mechanism. Estimated growth times and viscous damping rates for 
these unstable modes suggest that they are likely to lead to weaker instabilities 
than that of the fundamental r-mode !' This can be understood from the fact 



Uln this context, we should also mention that instabilities associated with the g-modes in 
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that a mode with a significant spheroidal part will be affected by (for example) 
bulk viscosity to a larger extent. 

A particularly interesting question, that provided the original motivation for 
the work of Lockitch and Friedmar£3, concerns what happens to all the r-modes 
of a non-isentropic star as As 0, i.e. in the isentropic limit. Recall that a 
non-isentropic star has an infinite set of r-modes for each I and m ^ 0, while an 
isentropic star retains only a single such mode for I = m. The general consensus 
at the present time is that most of the non-isentropic r-modes will change their 
nature by acquiring a spheroidal component as As — *■ 0, i.e. they become hybrid 
modes. Some light on the detailed behaviour has recently been shed by Yoshida 
and Leellij. They suggest that the r-modes of a non-isentropic star fall into 
one of three categories. The first class are the fundamental I — m modes (that 
have eigenfunctions Tu ^ a}~^ and lead to the strongest instability). These 
modes are largely unaffected by changes in As ■ Similar results were reported by 
Andersson, Kokkotas and SchutzcS. The second class of modes are the various 
I = m overtones ( whos e eigenfunctions have nodes inside the star). According 
to Yoshida and Leellia, these modes approach hybrid modes as the spin-rate of 
the star increases to a level where the slow-rotation approximation is no longer 
appropriate. Finally, the I ^ m r-modes of a non-isentropic star have frequencies 
that vanish (in the rotating frame) as the spin-rate increases. Consequently, 
these modes may not be relevant in a rapidly rotating star. 

It turns out that the general behaviour of the modes in the isentropic limit 
is rather similar to what happens at fast rotation rates. Until recently, low- 
frequency modes had not been calculated for rapidly rotating stars. The situ- 
ation has changed with calculations for the Maclaurin spheroidsE£j (essentially 



reproducing and elucidating very early work by Bryan££3) as well as for rapidly 
rotating polytropesEZI. These calculations are complicated by the fact that the 
mode structure changes as the spin rate increases, and for rapidly rotating stars 
one must include a large number of multipoles in a decription of ±he mode eigen- 
functions. Interestingly, the polytrope results of Yoshida et alKO indicate that 
the slow rotation expansion (taken to order Vl^) provides a good representation 
of the r-modes also in rapidly spinning stars. This means that the current mod- 
els for the spin evolution of young neutron stars etcetera may, in fact, be better 
than expected. 



7.2. Slowly rotating stars in general relativity 

One important step towards improving our understanding of the r-mode in- 
stability corresponds to describing the modes in general relativity. The Newto- 
nian picture of stellar pulsation, that we introduced in sections 2.1-2.3, is readily 
generalised to the relativistic case. The main difference is that the gravitational 
waves generated by the modes now appear as perturbations of the spacetime 



non-isentropic stars have been discussed by Lai 
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metric. Consequently, the various modes are no longer "normal modes" , but 
must satisfy outgoing-wave boundary conditions at spatial infinity. This pro- 
vides an additional complication in general mode-calculations (where the eigen- 
frequencies are now complex), but it is not relevant for a slow-rotation study of 
low- frequency oscillations. As can be deduced from ( p5[ ) the gravitational- wave 
dissipation will not enter until at very high orders in the slow-rotation expan- 
sion. Hence, the r-mode frequencies remain real to leading order, and we are 
(in principle) dealing with a "standard" normal-mode problem. 

Considerable progress towards a,n un der stan ding of the relativistic r-mode 



problem has been made recentlyBli^'E^jyiHSlHZl. Particularly relevant is work 



by Lockitch, Andersson and Friedmanlif3 that puts the various facets of the 



problem in the appropriate context and provides the first actual calculation of 
the relativistic r-modes. The picture that emerges with this work is somewhat 
more complex than the Newtonian one, but it can be understood if we recall 
the difference between isentropic and non-isentropic Newtonian perturbations. 
First of all, one can prove that (apart from a set of stationary dipole modesllHj) 
rotating relativistic isentropic stars have no pure r-modes (modes whose limit 
for a spherical star is purely toroidal). This is in contrast with the isentropic 
Newtonian stars which retain a vestigial set of purely toroidal modes (the fun- 
damental I = m modes that lead to the strongest instability). Instead, the 
relativistic corrections to the Newtonian r-modes with I = m > 2 are both 
toroidal and spheroidal. Thus these modes become discrete hybrid modes of the 
corresponding relativistic models. So far these modes have been computed for 
slowly rotating isentropic stars to first post-Newtonian order (~ M/R). 

As in the Newtonian problem, non-isentropic stars are somewhat different. 
In the slow-motion approximation in which they have so far been studied, non- 
isentropic stars have, remarkably, a continuous spectrum. Kojima EHJ has shown 
that purely axial modes would be described by a single, second-order differential 
equation. The continuous spectrum is implied by the fact that the correspond- 
ing eigenvalue problem is singular (the coefficient of the highest derivative term 
of the equation vanishes at some value of the radial coordinate) . The mathe- 
matical argurncnts for the continuous spectrum have been elucidated by Beyer 



and KokkotasE£3. However, as the latter authors point out, it is not yet clear 
that the continuous spectrum is physically relevant. In addition to the possible 
continous spectrum one can show that discrete r-mode solutions also exi.s t. Such 
discrete mode-solutions to Kojima's equation were recently calculatedH3. These 
modes are the relativistic analogue to the Newtonian r-modes in non-isentropic 
stars. 

Even though much work remains before the r-modes and their instability 
are understood in general relativity, the recent results provide much useful in- 
formation. In particular, even though the nature of the equations for isentropic 
relativistic stars is rather different from the corresponding non-isentropic ones, 
the particular modes that would be the analogues of the vestigial I = m r-modes 
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of isentropic Newtonian stars are not too dissimilarli23. In fact, the calculated 
mode frequencies differ only at the level of a few percent. The toroidal com- 
ponents of the velocity field are also virtually identical (although one must 
remember that the relativistic isentropic mode also has a significant spheroidal 
component). 

7.3. The role of the magnetic field 

Another important facet of neutron star physics that we have not yet dis- 
cussed is the magnetic field. Intuitively one would expect the interplay between 
a large amplitude pulsation mode and the magnetic field to be interesting. The 
simple fact that this may lead to observable effects is strong motivation for try- 
ing to incorporate the magnetic field in our models. To do this is, however, far 
from trivial and at the present time all results in this direction must be consid- 
ered as rather uncertain. On the other hand, the same is true for the r-mode 
spin-down scenario (section 5.1) so why should we shy away from speculating 
on the role of the magnetic field? 

There have so far been three main discussions of magniitic fields in the con- 



text of the r-mode instability. In the first of these, Spruitti2l outlined a possible 
scenario where an accreting neutron star with a weak magnetic field is spun 
up to the point where the r-mode becomes unstable (see Figure |lj). As we 
discussed earlier, it is likely that the r-mode instability then leads to a thermo- 
gravitational runaway that heats the star up and makes the mode-amplitude 
grow. Spruit suggests that gravitational radiation reaction then induces differ- 
ential rotation in the star. If so, the interior magnetic field could wind up until 
it reaches a critical point and becomes unstable due to buoyancy forces. When 
this happens, the star sheds a considerable amount of electromagnetic energy, 
perhaps as a gamma-ray burst, which slows it down significantly. Left over after 
the turmoil is a slowly rotating neutroii star with an extreme magnetic field (of 



the order of 10^'* — 10^^ G): a magnetaiQ'EI. This idea is undoubtedly speculative 
but it raises some interesting questions, and may also capture some features of 
the r-mode runaway. Qualitatively, the most important issue regards whether 
the instability leads to differential rotation or not. 

All the models that we have discussed in the previous sections assume that 
the star remains spinning uniformly. That this will actually be the case is 
far from clear. An important counterexample is the evolution of an unstable 
bar-mode in a rotating ellipsoid. As was first shown by MilleiH, an unstable 
mode drives a Maclaurin spheroid towards a Dedekind ellipsoid (more or less a 
rotating american football). En route the evolution is thought to proceed via 
a sequence of differentially rotating Rieniann ellipsoids. By analogy one might 
therefore expect that the r-mode instability will generate differential rotation. 
To conclusively prove this will, however, not be easy since it involves modelling 
how radiation reaction affects various parts of the star. 
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The question whether an unstable r-mode can be prevented from growing by 
the magnetic field was recently discussed by RezzoUa, Lamb and Shapiroli^. In 
an interesting approach to the problem, they considered the leading order result 
for the r-mode velocity field (basically the time-derivative of our equations ( ^8[ ) 
and (^9[)) as evolution equations for the fiuid elements. By evolving this system 
in time one finds a differential drift, the magnitude of which depends on the 
latitude of the fluid element and the mode amplitude. RezzoUa et al study how 
this deduced differential drift affects the magnetic field of the star (assuming 
that the field is frozen into the fluid elements). This leads to several suggestions. 
First of all, it is estimated that the r-mode oscillations cannot prevail if the star 
has a magnetic field ^ 10^^ G or higher. In a star with a weaker magnetic 
field, B > 10^" G, the rotation rate must be above 0.35^1k in order for the 
mode to survive. This means that the instability could operate both in young 
neutron stars {B ~ 10^^ G) and recycled ones {B ~ 10® G) provided that they 
spin fast enough. Secondly, the differential drift due to the r-mode twists the 
magnetic field and, just like in Spruit's model, this affects the strength and 
nature of the field. RezzoUa et al predict that, given an initial configuration 
with a strong poloidal field and a much weaker azimuthal one, the magnetic 
field evolves in such a way that after the r-mode spin-down the azimuthal field 
has been strengthened to the point where it is 4-8 orders of magnitude stronger 
than the (virtually unchanged) poloidal field. This would suggest that the r- 
mode instability generates strong azimuthal magnetic fields in young pulsars. 
This is an intriguing possibility, and it may be worthwhile investigating whether 
it could be tested observationally. 

These are interesting suggestions, but they come with a disclaimer. The de- 
rived differential drift is not determined consistently within perturbation theory. 
RezzoUa et al. use first order (in both a and fl) results to deduce a nonlinear 
effect (of order a^). In absence of true nonlinear calculations it is not clear to 
what extent the predicted drift will be present. 

Several other issues regarding the possible role of the magnetic field are 



discussed in a recent study by Ho and LaiE^il. First of all, they discuss the effect 
of the standard magnetic braking (that dominates the spin-down in absence 
of the unstable r-modes). From the standard results one finds that magnetic 
braking proceeds on a timescale 

B^iiW " ^ 10«Mi.4i?uf ^il'^-s s , (133) 

for a polytropic stellar model. Here B12 = B/10^^ G. By comparing this to 
the timescale of r-mode growth ( |65| ) we see that magnetic braking will only be 
competitive for neutron stars with B ^ 10^^ G and higher. In this context. 
Ho and Lai make an interesting observation. They note that, in absence of 
radiation, the canonical angular momentum is conserved. This means that if we 
imagine a situation where the star is spun down (or up) by some other agent, 
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then the amphtude of the mode must change in order to keep Jc constant. 
Hence, one would expect to find a ^ This would not be important 

under normal circumstances, but it could be highly relevant if other spin-down 
torques dominate, eg. for a young magnetar. 

As we have already discussed in Section 3.1, a mode may be driven unstable 
by any radiative mechanism. Hence, it is interesting to try to quantify the extent 
to which electromagnetic radiation from the r-modes may affect the growth rate 
of the instability. Ho and Lai estimate this effect, and based on their results we 
find 

torn ~ -4.8 X 10l"7\fi.4Bi2^i?ro^Pl3 S . (134) 

Clearly, this effect would be relevant only for very strong magnetic fields. Basi- 
cally, the electromagnetic driving of the mode would become competitive with 
gravitational radiation for B 10^^ G, so again we deduce that the effect could 
be of importance for magnetars. Finally, Ho and Lai also point out that the 
shaking of the magnetic field lines due to the mode-oscillations will generate 
Alfven waves in the magnetosphere. These can also drive the mode, and one 
can estimate that 

tA ~ -2 X 10^ Ml, iR^oB-^^ s . (135) 

Interestingly, this leads to a more significant driving of an unstable mode than 
direct electromagnetic radiation (134). 

The above results suggest that the magnetic field may have great effect on 
the r-modes in magnetars, i.e. for stars with magnetic fields above (say) 10^^ G. 
For stars with weaker magnetic fields, such as recycled neutron stars in LMXBs 
the picture is not so clear. The outcome depends crucially on whether the r- 
mode leads to differential rotation/drift in the fluid or not. There is also the 
obvious caveat regarding neutron stars with a solid crust. It could be that the 
magnetic field is essentially frozen into the crust and that waves in the fluid 
core do not significantly penetrate the crust. If this is the case, the magnetic 
field may not have a large effect on the r-modes at all. On the other hand, it is 
expected that the core fluid will be electromagnetically locked into corotation 
with the crust on a very short time-scale. This may then prevent oscillations 
in the core entirely. Much further work into the interplay between neutron star 
pulsation and the magnetic field is needed before we can claim to have even a 
rudimentary understanding of this problem. 



7.4. Nonlinear effects 

It is clear that most of the results we have discussed in this article may 
be strongly affected by nonlinear effects. Nonlinear stellar pulsation theory is a 
relatively unexplored research area that provides many conceptual and technical 
challenges. Still, if we want to understand the actual evolution of an unstable 
pulsation mode and possible effects it may have on the bulk rotation of a star 
we must explore the non-linear regime. 
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In the case of the r-mode scenarios discussed in Section 5, one main unknown 
parameter is tlie saturation ampUtude as . The values we used in our discussion 
(as = 0.01 — 1) were chosen ad hoc, by assuming that the r-mode would be 
able to grow to a "large" amplitude. But it is not at all clear that this will be 
possible. Several nonlinear effects may set in at considerable smaller values of a 
and prevent further mode-growth. If this is the case, the various astrophysical 
scenarios must be altered accordingly. Some candidates for nonlinear mode 
saturation are: i) coupling to other pulsation modes, that could lead to a cascade 
of energy from the unstable mode into other modes, ii) turbulence, which is likely 
to play a role since the Reynolds number is very high in the r-mode problem, 
and iii) nonlinear effects that change the nature of the mode in such a way that, 
for example, shocks develop. At the present time, the only available results 
regarding the nonlinear mode saturation follow from an eaiixnate of turbulence 



that arises in the Ekman layer at the base of the crustlll3. These estimates 
suggest that the shear motion in the Ekman layer becomes turbulent when the 
mode reaches an amplitude** 

ac ~ 1.6 X lO-^A-^Tg^^P%^ (136) 

where A — |Ai;|/|i;| as in Section 4.5. We recall that a typical value might beS 
A « 0.1. Hence turbulence would play a role provided that the mode can grow 
to an amplitude a « 0.1. Since the turbulent dissipation rate increases as the 
cube of the amplitude while the gravitational-radiation growth scales as a^ one 
can estimate the the mode will saturate whentiJ 

as « 3.5 X 10-3A-3pr| (^^^p) (137) 

where Cd is the relevant drag coefficient. From this result we can deduce that 
turbulence may saturate the r-mode growth at amplitudes significantly below 
unity provided that A is relatively large. But if A « 0.1 we get as « 3.5 and it 
seems likely that other nonlinear mechanisms will provide the mode-saturation. 

A second issue of utmost importance regards the spin-evolution of the star 
during the instability phase, i.e. the back-reaction of the mode onto the bulk 
rotation of the fluid. In the phenomenological spin-evolution model discussed 
in Section 5.1 it was assumed that, once the r-mode had reached the saturation 
amplitude, all the excess angular momentum was drained from the bulk rotation. 
That this will actually be the case is far from clear, and would in fact seem 
rather unlikely. Nonlinear effects may well cause differential rotation in the 



*It should be noted that the turbulence estimates are based on experimental data for the 
flow of water above rugged surfaces. It is not clear to what extent such results will be 
relevant for neutron star oscillations. One should also note that it is assumed that one can 
meaningfully use the standard r-mode gravitational wave estimates even in the turbulent case. 
This corresponds to assuming that the turbulence is confined to the Ekman layer, which may 
be dubious since turbulence is by its very nature convective. 
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star (or a secular drift as in the study of RezzoUa et alll^), and if this is the 
case our current modeUing will be inadequate in many ways. A first attempt 
to gain insight into the non-linear back-reaction has been made by Levin and 



UshomirskyEid. The considered the toy-problem of an r-mode in a shell, and 
find that the r-modes do indeed generate differential rotation in the star. 

It seems to us that we need numerical studies using fully nonlinear hydrody- 
namics to investigate these difficult issues. This is obviously a very challenging 
task. Ideally, one would want to study the onset of instability and follow the 
mode through to saturation and then see what effect the instability has on the 
spin of the star. Computationally, this means that one would need a numer- 
ical evolution that resolves the mode-oscillation (on a millisecond timescale) 
and tracks the star through hundreds of seconds. This is certainly not possible 
given current technology. As we understand it, the best one can hope for at 
present is to study the star over (say) a few tens of rotation periods. This will 
not allow a complete study of the problem, but it may still provide answers to 
many crucial questions. Efforts in this direction are underway, based on either 
direct numerical evolution of the genera l-relativistic hydrodynamical equations 
on a fixed background spacetimct23lH3 or as an integration of the Newtonian 
hydrodynamical equations and the use of an " accelerated" gravitational radia- 
tion reaction forcelHj. As we were finishing this review the first results of these 
efforts became available. Simulations due to Stergioulas and Font indicate that 
the r-modes may be able to grow to a surprisingly large amplitude. In fact, 
there are no significant signs of mode-saturation until at unrealistically large 
amplitudes. These results are very promising, and it seems likely that contin- 
ued efforts in this direction will provide important insights into the nonlinear 
physics of the r-mode instability. 



8. Concluding remarks 

In this review we have discussed the r-modes in rotating neutron stars and 
the associated gravitational-wave driven instability. This is a research area 
that has attracted considerable interest following the prediction that hot young 
neutron stars may be spun down to a level comparable to that extrapolated 
from observations and the suggestion that the associated gravitational waves 
may be detectable with the new generation of interferometers. 

Our aim in writing this review was to summarize most of the ideas and sug- 
gestions that have been made regarding the r-mode instability and put them 
in the appropriate context. In doing this we obviously had to make sacrifices, 
but hopefully we have managed to provide a useful introduction to the many 
relevant issues. What should be clear is that our present understanding of in- 
stabilities in rotating neutron stars and their potential astrophysical relevance 
is unsatisfactory in many respects. All current suggestions are based on simpli- 
fied, often phenomenological, models. The hope must be that future investment 
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in this field will lead to a truly quantitative description of "realistic" neutron 
stars, and that the various proposed astrophysical scenarios will be described 
in detail rather than "order of magnitude" . We have far to go before we reach 
this goal, and it may well be that actual observations will beat us theorists to 
the answers of the many challenging questions that we are just beginning to 
formulate. 
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